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PREFACE 


It is widely felt that when the equations of mathematical 
physics are written out in full Cartesian form the structural 
simplicity of the formulae is often hidden by the mechanical 
labour of writing out every term explicitly. Attempts have 
been made to reduce this labour by one form or another of 
vectOT algebra; but it has always seemed to me that this 
method both introduces new difficulties and is insufficiently 
general. Thus the product of two vectors, in vector lan- 
guage, means one of two things, either the scalar or the 
vector product, and it is not physically obvious why just 
these functions of the vectors should arise and no others. 

The use of tensor notation, with the summation con- 
venticm, carries out as great a simplification of the writing 
as does vector notation. The notation has actually at- 
tracted attention owing to its applications in the theory of 
relativity, but for ordinary purposes two great abbrevia- 
tions may be made. We use rectangular Cartesian axes; 
the result is that the distinction between covariant and 


contravariant vectors disappears, and with it the terms 
arising from curvature of the surfaces of reference. The 
formidable character of most of the formulae of the theory 
of relativity is absent from the formulae of tensors referred 
to Cartesian axes. The tensor method is a necessity for re- 
lativity; for applications in dynamics, electricity, elas- 
ticity, and hydrodynamics it is a great convenience. 

It is found that the scalar and vector products are not 
the only functions of two vectors that arise, though the 
theory provides reasons why they are important in many 
applications. There is also a symmetrical product, which 
ordinary vector notation is completely unable to express. 
In tensor notation it arises naturally as a symmetrical 


VI 


PREFACE 


tensor of the second order. The system of moments and 
products of inertia of a rigid body constitutes such a 
tensor; so do the stress components and the strain com- 
ponents in an elastic solid. 

The present method, like vector notation, is of use prin- 
cipally in proving general theorems. In concrete applica- 
tions there is usually some asymmetry about the coor- 
dinates that makes it necessary to abandon the tensor 
form at some stage in the work. It has been said that 
vector equations are like a pocket map, and it has been re- 
plied that a pocket map has to be taken out of the pocket 
and unfolded before it is of any use. The same applies to the 
tensor method, and for the same reason; but it has the 
great advantage that it is not a new notation, but a concise 
way of writing the ordinary notation, so that the unfolding 
can be carried out more conveniently when occasion arises. 

What is usually called Statics is treated in Chapter v, 
after Dynamics. I consider this to be the proper order, be- 
cause Statics is a special case of Dynamics, and many of its 
formulae have physical significance for reasons explained 
in Dynamics. The customary reversal of the order is due, 
I believe, to the fact that an introduction to mechanics has 
to be given at schools before the students have received 
any training in calculus; but this need not influence 
students working for a university examination. 

It should perhaps be stated that the object of this work 
is to illustrate the use of tensor methods; it does not claim 
to give a complete theory of all the subjects touched, re- 
ference for which must be made to the standard text-books. 

I must express my gratitude to ]\Ir M. H. A. Newman, 
Miss L, M. Swain, Dr S. Goldstein, and Dr Bertha Swirles 
for assistance at various stages in the work, and to the 
staff of the University Press for their care in the printing. 

HAROLD JEFFREYS 


September 1931 



NOTE 


Since this book was written most of the material in it 
has been incorporated in Chapters 2 and 3 of Methods of 
Mathematical Physics^ by my wife, Bertha Swirles Jeffreys, 
and myself. The present reprint has been made partly 
because of a continuing demand for a treatment of Car- 
tesian tensors by themselves ; partly because some results, 
notably on the thermodynamics of an elastic solid and the 
circulation in viscous flow, are not given in textbooks of 
the special subjects. 
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In the present reprint a few small changes have been made. 
In Ex. 3, p. 15 j the components are given explicitly; these 
have been found useful in developing the strain energy for 
an elastic sphere under rotation. Omissions have been 
corrected on pp. 81 and 82. As these do not affect the 
later work they appear to have given readers no trouble. 
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December, 1968 


CONTENTS 


Chap. I. Cartesian Tensors 1 

II. Geometrical Applications 16 

III. Particle D 3 mamics 24 

IV. Dynamics of Rigid Bodies 29 

V. Equivalence of Systems of Forces 42 

VI. Continuous Systems 48 

VII. Isotropic Tensors 66 

VIII. Elasticity 71 

IX. Hydrodynamics 83 

Index 93 



If we have two sets of rectangular axes (Ox, Oy^ Oz)y 
{0x\ 0y\ Oz*) at the same origin, the coordinates of a 
point P with respect to the second set are given in terms 

of the coordinates with respect to the first set by the 
equations 


x' = l^x + rri^y + n^z 
y* ~ l^x + m^y + n^z 

2' = ZaX + m^y + n^z 



The quantities (^, Wi, , ng) are the cosines of the 
angles between the various axes ; thus is the cosine of the 
angle between the axes Ox' and Ox, is the cosine of the 
angle between Oy' and Oz, and so on. It follows that the 
coordinates {x, y, z) can be expressed in terms of (a:', y\ z') 
by the relations 


X — Ij^x' + l^y’ + l^^z* ^ 

y = niyx' + m^y' + m^z ' t ( 2 ), 

z = n^x' + n^y' + ^ 32 ' j 

We can shorten the writing of ( 1 ) and ( 2 ) considerably 
by a change of notation. Instead of (x, y, z) let us write 
(^iX^y ajg), and instead of (a:', y\ z') write (x/, x^y x^). We 
can now say that the coordinates with respect to the first 
set of axes are x^y where i may be 1, 2, or 3; and those with 
respect to the second set are x^ y where may be 1 , 2 , or 3 . 
Then in ( 1 ) each coordinate x^' is expressed as the sum of 
three terms depending on the three Xi . Each x^ is associated 
with the cosine of the angle between the direction of that 
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Xi increasing and that of x/ increasing. Let us denote this 
cosine by Then we have, for all values otj, 

Xi = a^iXi^ + a^iX^ + a^iX^ 

= S dijXi (^)* 

t«l, 2, 3 

Conversely (2) can be written 

S a^x/ (4), 

y=i.2. 3 

the a^i having the same value as in (3), for the same values 
of i and 7 , because it is in both cases the cosine of the angle 
between the directions of Xi and x/ increasing. 

In mathematical physics we often have to deal with sets 
of three quantities in relation to a set of axes, of the general 
form Ui (that is, tq, u^y v^)y and such that in relation to a 
different set of axes the corresponding quantities are 
(Wj', u^y u^')y which satisfy the relations 



Uf ~~ L a^fUi 

i“1.2, 3 

(5) 

and 

Ui= S CLifU/ 

(6). 


3 = 1 , 2 . 3 

Such sets of three quantities are called tensors of the first 
order y or vectors. The individual v^yU^yU^ may be called the 
components of the tensors. 

Clearly if we multiply all of the Uf and u/ by the same 
quantity m we get 

mu/ = S dii (mUi) (7), 

i-l. 2, 3 

so that mUi is another tensor of the first order. 

Again, if we have two tensors of the first order, w,- and 
Viy we shall have 

u/ + v/ = S O..^ (m,. + Vi) (8), 

<•=1, 2. 3 

SO that Uf + is a tensor of the first order. 

We notice that each of the equations (3) to (8) is really 
a set of three equations ; where the suffix i or j appears on 
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the left it is to be given in turn all the values 1, 2, 3, and the 
resulting equation is asserted in each case. In each such 
equation the right side is the sum of three terms, obtained 
by giving^ or i the values 1, 2, 3 in turn and adding. Wher- 
ever such a summation occurs a suffix is repeated in the 
expression for the general term; where there is a summa- 
tion for all values of the general term, such as a^u/y con- 
tains^' twice. We make it a regular convention that, unless 
the contrary is stated, whenever a suffix is repeated it is 
to be given all possible values and that the terms are to be 
added for all. Thus we write (5) as simply 



u/ = atfUi (9), 

the summation sign being automatically understood by 
our convention. Then (9) really means three equations, 
with three terms on the right of each, but we can by means 
of our conventions express all of the twelve terms com- 
pactly by the single equation (9). 

There are single quantities, such as mass and distance, 
that are the same for all sets of axes. These are called 
tensors of zero order y or scalars. 

Consider now two tensors of the first order, Ui and v*. 
(When we write “a tensor we mean of course a tensor 
of the first order whose components are This is 

another piece of shorthand.) Suppose each component of 
the one multiplied by each component of the other; then 
we obtain a set of nine quantities expressed by where 
oach of i and k is independently given all the values 1, 2, 3. 
The components of Uf , v* with respect to the other set of 
es are u/, Vi say ; and 

u/vi' = (atfUi) 


= afia^iUiV^ ( 10 ). 

A he suffixes i and k are repeated on the right. Thus (10] 
represents nine equations, each with nine terms on the 
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right. Each term on the right is the product of two factors, 
one of the form depending only on the orientation of 

the axes, and the other of the form representing the 
products of the components referred to the original axes. 
In this way the various u/ v{ can be obtained in terms of 
the original But products of two vectors are far from 

being the only quantities satisfying this rule. In general a 
set of nine quantities Wij^ referred to a set of axes, and 
transformed to another set by the rule 

( 11 ), 

is called a tensor of the second order. 

We may go on similarly to construct and define tensors 
of the third, fourth, and higher orders. Thus a set of quan- 
tities that transforms like is called a tensor of 

order n, where n is the number of factors in this product. 

When we say that a certain set of quantities is a tensor 
of any order n, we mean that we have ways of specifying 
its components with respect to any set of axes, and that the 
components with regard to any two different sets of axes 
are related according to the rule appropriate to tensors of 
that order, and in particular to the products of the coor- 
dinates with n factors. For instance, if we say that is a 
tensor of order 1, we are not simply defining u/ as meaning 
We are supposing both that uf has a meaning, such 
as a displacement or a velocity, with reference to the axes 
of xf, and that the value of each component is equal to 
aiiUi. Thus the statement that any set of quantities is a 
tensor is not a mere convention, but a statement capable 
of test and therefore needing proof. In (7) and (8), for 
example, our data are that w,- and u/ are the components 
of a vector with regard to two different sets of axes. We 
prove that the sets of quantities obtained by multiplying 
both by the same quantity are related according to the 
vector rule; and therefore the products are vectors. 
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If we interchange^' and Z in (11), we get 

( 12 ). 

But on the right % and h are “dummy suffixes”; that is, 
they are to be given all possible values and the results 
added. It is ummportant which of them we call % and which 
h\ we may therefore interchange them and get 

Wi! = akiatjWf,( = ( 13 ). 

Thus transforms according to the same rule as and 
therefore is another tensor of the second order. The im- 
portance of this is that if we know the set of quantities 
arranged 

/tVii 

( ^21 ^22 ^23 ) ( 14 ) 

\uj 31 MJ 32 yj^/ 

to be a tensor of the second order, tlien the arrangement 

/ W^31 \ 

I y>X2 ^22 1^32 I (15) 

' ^^23 ^33 

is another tensor of the second order. Therefore the sets 
{y^ik + y^ki) and (Wi^ ~ Wfit) are tensors of the second order. 
The first of these has the property that it is unaltered by 
interchanging i and k, and is therefore called a symmetrical 
tensor. The second has all its components reversed in sign 
when i and k are interchanged, and is called an antisym- 
metrical tensor. Clearly in an antisymmetrical tensor the 
“leading diagonal” components, i.e. those with i and k 
equal, are all zero. Also, since 

^ik = i (Wif, + Wkt) + i (w^ifc - ^ki) ( 16 ), 

we can consider any tensor of the second order as the sum 
of symmetrical and antisymmetrical parts. 

The gradient of a scalar is a vector. For if U is a scalar, 
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its gradient is dUjdXi or dUjdx/ according to the set of 
axes. But 


dU dxi dU dU 

dx/ ~ dx/ dXi dXi 


(17), 


so that the gradients transform according to the vector 
rule. Similarly the gradient of a vector is a tensor of 
order 2. For if Uf and u/ are the components of a vector 
with respect to two sets of axes, 


du/ _ dxfi du/ 
dxi ” dxi dx^ 



dui 


( 18 ), 


so that the rule of transformation is as in (11). 

Since Xf is a vector, it follows that dXijdXk is a tensor of 
the second order. But dxfjdxf, is unity ii i = k and zero if 
i ^ k. Hence the set of quantities 8,*, such that 



constitutes a tensor of the second order. We can prove this 
directly; for if we apply (H), in the new system of 
coordinates should be given by 

8// = (19). 

The suffix k has to take all values 1, 2, 3. But M k ^ 
8rt is 0, and the corresponding term is zero. If A: = i, 
S,fc = 1, and the result of the summation with regard to A: is 

8 ^/ = (lijUii ( 20 ). 

But the are the direction cosines of the axis of x/ with 
regard to the Xi^ and the an are those of x/ with regard to 
X, . Hence a^an is the cosine of the angle between x/ and 
x/, and is equal to 1 if the axes are identical and to 0 if they 
are perpendicular. It follows that the result of the trans- 
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formation is that 8,,' = 1 if J = ;, and 8,,' = 0 if Z. It 
foUows tli&t the set of titles 


/I 0 0\ 

I 0 1 0 ) (21) 

\0 0 1 / 

is transformed into itself by the rule ( 1 1 ) and therefore is 
a tensor of the second order. 

n is a vector and we form the product 8,„u ^ , we have 
a tensor of order 3. But now put m Jk and add for all 
values of L Since S,* = 0 except for ^ = t, the only term 
different from zero is that for k = i, and this is u,. Hence 


^ Ui ( 22 ). 

This operation therefore replaces the suffix k by i. The 
tensor can therefore be called the substitution tensor. 

In the tensor m;,,. let us put k = i, and in accordance 
with our convention add for aU values of i. Then the 


corresponding quantity 
summing; but 


is got by putting and 


Wf/ — 

= Wii (23). 

Thus transforms into itself and therefore is a scalar. 

This operation of putting two suffixes in a tensor equal 
and adding accordingly is known as contraction. In general 
it gives a new tensor, whose order is less by 2 than that of 
the original tensor. If for instance we contract the tensor 




we obtain 


UiVi *= + U2V2 + U2V2 




which is the scalar product of Ui and . 

Similarly the tensor yields the symmetrical and 

antisymmetrical tensors (u{v^ + Uj^Vi) and (u-Vf^ — u,^Vi). We 

may call these the symmetrical and antisymmetrical products 

of Ui and 
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The tensor dujdxi gives similarly, on contraction, a scalar 

dUi 
dxi 


dill du2 dit^ 
dxo dx. 


dXi ^*^3 


(25), 


which is known as the divergence of Ui ; while it gives also 
symmetrical and antisymmetrical tensors 

dui 


dXi dxv 


and 


dXi 


The former has important applications, especially in the 
theory of elasticity and hydrodynamics; the latter is 
known as the curl or rotation of The vanishing of the 
curl is the condition that w, may be the gradient of a scalar. 

All the above considerations can be extended to any 
number of dimensions. In n dimensions a tensor of order r 
has components. A tensor of order 2, in particular, has 
components. If it is antisymmetrical, the n diagonal 
components are zero, and the others are equal and opposite 
in pairs. Hence an antisymmetrical tensor of order 2 has 
^n{n — 1) independent components. It n — 1, 2, 3, 4,... in 
turn, this number is 0, 1, 3, 6,... . It happens that in three 
dimensions the number of numerically independent com- 
ponents of an antisymmetrical tensor of the second order is 
equal to the number of components of a vector. Actually it 
can be proved that with any vector we can associate an anti- 
syrametrical tensor of the second order, and conversely. 
This is not true in any number of dimensions other than 3. 

Since the a.^ are the direction cosines with respect to the Xi 
of three perpendicular lines, they are connected by six re- 
lations 



-f- ^21^ 

+ = 

r 



ai2^ + 022^ 

+ 032* = 


(26), 


"b ^23^ 

+ 033*' = 

i] 



-f- Cf22^23 


^32^33 “ 

0 1 



-f 023^21 

+ 

^33^31 “ 

0 

(27). 

^11^12 


+ 

^31^32 ~ 

0^ 
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We notice that the second and third of (27) both contain 

(®ii> « 3 i)* We may therefore solve them for the ratios 

of these quantities. Thus 


Oil 


®2I 


a 


31 


®33®22 ^®12 “ ®13®32 ®13®22 ~ ^23®12 


- k 
(28), 


say. Substituting in the first of (26) we get 

1 = A: {Oii ((X23®32 ~ ®33®22) + ®21 (®33®12 ~ ®13%2) 

“i" ®31 (®I3®22 “ ®23®12)} 


= - A; 


Also 


®21 

^2 ®22 ®32 
®13 ®23 ®33 


(29). 


^^{(®23®32 ^^22)^+ (®33®12~ ®13®32)^ + (®13®22 ~ ®23^12)^} = i 

^ (30). 

But we have a general identity 

(a^ + 62 + c2) (a'2 + 6'2 + c'2) - (oa' + 66' + cc')2 

= (6c' — c6')2 + (ca' — ac')2 + (06' — 6a')2 (31), 

Hence 


k^ [{aj- + 4 . 

(®12®13 4" ^22^^23 ^ 32 ^ 33 )^] ~ 1 (32). 

But on account of the second and third of (26) and the first 
of (27) the expression in brackets is unity, and therefore 

k = ±l (33). 

For any given transformation the determinant in (29) is 
therefore equal to ± 1 . Evidently its sign is reversed if we 
interchange any two of the suffixes j, for this interchanges 
two rows of the determinant; so that the sign is a matter of 
the numbering of the axes. If we start with a rigid frame 
attached to the axes a:,, and rotate it continuously till it 
is attached to the axes x/^ all the vary continuously and 
therefore the determinant cannot change from + 1 to — 1 

a r^T 
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or from — 1 to + 1. If then goes to to x^y and Zg 

to x^y the determinant is initially 

1 0 0 = 1 (34), 

0 1 0 

0 0 1 

and therefore 

A: = - 1 (35), 

and the determinant formed by the is always + 1. 

If we have a frame of axes {x^x^x^^ we can turn it by a 
continuous movement so as to bring x-^ along the old x^y 
x^ along the old x ^ , and x^ along the old Xi . In this case we 
have 

Xj = X 2 ~ X^y X^ ~ Xi (3®)i 

and 

CEu = 0, ~ 1 > ®31 “ ^ > ®12 “ ®22 ~ ^32 “ f » 

®13 = 1» ^23 “ ®33 ~ ^ 

The determinant of the a,y is therefore 

0 10=1 (38), 

0 0 1 

1 0 0 

as before. Any rotation of the axes that does not alter the 
cyclic interchange of suffixes 1 2 3 1 2 3 1 ... therefore 
leaves the determinant equal to unity, and therefore so 
long as we always use right-handed or always left-handed 
axes the determinant of the a^y is H- 1. 

With this restriction 



flii — ^22^33 


32^23 > 


Qoi ®13^32 


33^12 > 


®31 — ^ 23®12 ® 13 ®; 


(39), 


and therefore every direction cosine is equal to its first 
minor in the determinant. 
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These relations are of course identical with those ex- 
pressed in the usual notation of solid geometry by 

(40). 

Now suppose that is a vector, and consider the set of 
quantities 

/ 0 «3 - ^\ ( 41 ). 

— 0 

U 2 — 0 

■^PP^y (11) to this, taking i to be the number of the row and 
k that of the column. We see that Uj enters as and as 
— Its coefficient in Wn* is therefore 
and in all 

Wji = ®2l) 

+ (fli/Uai — Cl2i^u) ^ (42). 

This is obviously zero iij — 1. If Z and if the other axis 
perpendicular to x/ and x/ is x^ , oxid jlnjln is a cyclic 
order, the quantities in brackets are equal to (a^, a 2 „, a^n)- 
This is true if I immediately succeeds j in the order. If I 
precedes j by one place the signs are reversed. Hence if 
j = 1 and Z — 2, or j = 2 and Z = 3, or if ^* = 3 and Z = 1, 


= CLin^i = w„' 

and in the alternative case 


(43), 


Thus 





= — u 


u. 


n 


— U 



(44) . 

(45) , 


0 

Uz — Ui 

and is of the same form as (41). Thus with any vector we 
can associate an antisymmetrical tensor of the second order. 
Conversely with any antisymmetrical tensor of the second 
order we can associate a vector. 
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We can proceed alternatively by considering the set of 
quantities defined by the condition that if any two 
of iy hy m are equal the corresponding component is 0 ; if 
iy ky m are all imequal and in cyclic order, the component 
is + 1 ; if the order is not cyclic, the component is — 1 . 
Let us see whether this is a tensor of the third order. If so, 
we should have 

^iln = 

(^®)* 

Now if, for instance,^* = I, the right side is clearly zero and 
^iin = 0. If jy ly n are all unequal, the expression is 







(47), 


®2n 

®3n 


which is equal to 1 if are in cyclic order and to — 1 if 
not. Hence the set of quantities is transformed into 
itself by the rule for transforming tensors of order 3, and 
therefore constitutes a tensor of order 3. This is called the 
alternating tensor. 

Now consider the product where Uj, is a vector, 

Tliis is a tensor of the fourth order. If we contract it by 
putting p = m and summing we get a second order tensor 
^ik = • If i = 1 and k ~ 2y the only value of m that 

makes €i 2 mdifferentfromzerois 3 ,and then + 1. Hence 

«Ji 2 = W 3 (48). 

If t = 2 and k = 1 , w is 3; but 213 is the reverse of cyclic 
order and €213 = — 1. Hence 

Similarly we find that the elements of are 

/ 0 «3 - W2\ 

( - «3 0 uA 

\ 2^2 "■ / 


( 50 ), 
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so that the antis 3 mimetrical tensor associated with a vector 
can actually be obtained from it by multiplying by 
and contracting. 

Again, suppose that we are given a tensor of the second 
order and that we form a vector by multiplying by 
and contracting twice. We have, if m = 3, 

^ — ^ ikS^ik ~ ^123^12 ^213^21 ~ ^12 *^21' 

Thus — Cijcjn'^ik (^f)* 

If Wijf is symmetrical, this evidently gives zero. If it is 
antisymmetrical the components of are numerically 
twice those of . 

On accoxmt of the intimate relation between the vector 
and the antisymmetrical tensor we shall habitually denote 
the tensor of (41) by so that 

^11 ~ ^22 “ ~ ^ * "^2 ~ ^ * ^23 ~ "^31 “ ^2 ’> 

^21 “ ^ I ^2 “ » ^3 “ ^2 

It will always be seen at once whether the vector or the 
tensor is intended, since the former has one and the latter 
two suffixes. 

If we have any three vectors , w;» , and consider the 
scalar we see that 

— U^ViW^ — U-^V^tOi — UiVqW^ 

W3 (53), 

Vl ^2 

so that we have a concise way of writing the determinant 
formed by the components of three vectors. If any two of 
the vectors are parallel this scalar vanishes. 

In associating a vector with an antisymmetrical tensor 
of order 2 a sign convention clearly arises. We make the 
positive signs in (41) lie one place to the right of the leading 
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diagonal. If then we have two vectors and v*, their anti- 

symmetrical product is u^v,, — Uj,Vi, and in the associated 

vector we give the positive sign to when k follows i in 

the cyclic order. Hence the components of this vector are 
taken to be 

This vector is perpendicular to both the original vectors ; for 

Ui (u^V3 - V^V^) + Wj (ligUj — v^v^) -I- ^ 0^ 

('>hVa - (u^Vi - tC^Vj) + V 3 {WyV^ - U^Vi) = 0. 

We caU it the vector product of m,. and , and can save 
writing by denoting it by [u, . 

Similarly with the antisymmetrical tensor ^ we 

dXi dXk 

associate a vector so as to leave the sign unaltered when 
k follows i in the cychc order. Thus the components are 

du^ du^ 9w2\ 

Wj 3x3 ’ 9 x 3 ~ax,’ 0xi~a^/ 

This is often called “curl u.^* 

The Tensor Since this tensor is the product of 

two third order tensors, once summed, it is a tensor of 

the fourth order, i, k, m, p being arbitrarily assignable. 

Evidently \i i ~ k or m = p^ the corresponding component 
is zero. 

If i = m, the contribution from any value of s is zero 
unless also k = p^ and then 

^mp$ ~ i 

and the component is + !• 

If i = p, then no value of 5 gives a contribution unless 

k=m. Then one of e,*, and is + 1 and the other - 1 , 

and the component is - 1 . Hence the components of the 
tensor are as follows. 
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Ji i ^ rriyk = 'Py the component is + 1, unless i also = ky 
% — Pyk =my the component is — 1, unless i also = k, 
i k OT m = py the component is 0. 

These results apply also to the tensor 


and therefore 




^iks^mps ^ip^km 


(55). 


EXAMPLES 

1. If u/y Ujg" are the components of a vector with regard to 
three sets of axes, prove that the values of Uj^" are the same as 
would be obtained by transforming first from to u/ and then 
from u/ to Ujf'. 

2. Prove that Sj-,- = 3; “= 0. 

3. Show that 

+ Sa„8j,-Scjfc Sam^bk^ci’ 

4. Prove that 

Prove that ^tka^mps — ^aik^amp “ 

6. Prove that if u^, Vjg, are vectors, 

[w [w, v]]^ = (u.w,) - v„ iu^Wf), 


7. If 

A (u) = 

«ii 

«12 

^18 




W22 

^23 



«31 

^32 

^33 


prove that 


^ifk A (w) = ^imn^a^Jm'^kn* 

*iik‘^n'^im'^kn “ ^Imn A (^)» 

6A (w) = 

8. Use Ex. 7 to prove the rule for the multiplication of deter- 
minants 

A («) A (u) = 6A (u) A (v) = 6A {uv), 
where (wv),p = 


C H A FT E R II 


GEOMETRICAL APPLICATIONS 

The displacement from any point to any other obviously 
constitutes a vector. The distance between the points is a 
scalar. If Xi , y,- are the coordinates of the points and r the 
distance between them, 

r2 = (yt - Xf)^ (1), 

the square on the right indicating the scalar product of the 
vector into itself. Also the quantities (y^ — ar^/r constitute 
a vector. 

If we take a fixed point and consider points given by 

= a* + hr (2), 

where r is a variable scalar and the are constants such 
that 

= 1 ( 3 ), 

(X, - = r* (4), 

so that r is the distance of Xi from . If we take another 
point yi such that 

+ h^ (5), 

(Vi - (6), 

(Vi - = (-5 - r)2 (7), 

and therefore the distances between a,-, a;,-, and y^ are such 
that the sum of two of them is equal to the third. Thus the 
points are on a straight line; and (2) gives the equations 
of the line in terms of the parameter r. The are the 
direction cosines of the line. 

If we take two lines through a, given by 

Xi = ai-\- hr 

yi = ai + JUiS 


( 8 ). 

(9), 
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GEOMETRICAL APPLICATIONS 
the distance between Xf and is given by 

(y< - = {snii - rli)^ 

= 52 _|_ ya _ 2rslimi (10). 

But this quantity is also equal to ~ 2rs cos 6, where 

6 is the angle between the lines. Hence the angle between 
two intersecting lines is given by 

cos 0=liin( (11). 

If two lines have the same direction cosines they are said 
to be parallel. If two lines do not intersect we can take a 
line through any point on one of them parallel to the other ; 
then this line is inclined to the first at an angle given by ( 1 1 ) . 
We can then use (11) to determine a unique quantity 
associated with any two lines, which we may call their 
inclination, whether they intersect or not. 

If we have a line given by (2) and yi is a point outside it, 
the line joining and subtends a right angle at Xi if 

(Vi - ai)^ = (Xi - ai)^ + (Vi - XtY 

= + (y, - a, - rl,)^ 

= (Vi - %)* - (Vi — (12), 

and therefore 

^ = li {Vi ~ ^i) (13)- 

This gives the projection of the displacement y< — on 
the line. The foot of the perpendicular is 

+ hr = a, + hl^ iVk - «*)• 

Evidently r in (13) will be the same for all points y^ such 
that hVi is constant. Hence 

hVi = ^ (14) 

represents a plane perpendicular to the line. 

If we take two intersecting lines given by (8) and (9), we 
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can find the equation of the plane containing them as 
follows. If this plane is 

n^Zi = p (15), 

this equation must be satisfied by Xi and for all values of 
r and a. Hence 

riiat = p (16), 

riih - 0 (17), 

niTTii = 0 (18), 

and from (15) and (16), 

n, {Zi - ad = 0 (19). 

Then (17), (18), (19) are three homogeneous equations in 
the riiy and can be consistent only if 

i^i - «,) = 0 (20). 

This is the equation of the required plane. Also the Uf are 
proportional to 

~ (^ 2 ^^ ^3^2 > ^3^1 — ^ 2 ^ 1 ) (^1)* 

But n..2 = 1 (22), 

(/gma - ^3^2)2 + (hrrvi - hrn^Y + (km^ - 
= + l^) + m^) — (l^m^ + 

= 1 — cos^ 6 

= sin2 e (23). 

Thus sin = ± 7re„. (24). 

The ambiguity in sign corresponds to a general one in 
specifying the parameter r of a point on a line. If r in (8) 
is taken negative, we get a point on the line on the opposite 
side of ai from those given by positive values of r. But if 
we reverse both r and the li we still keep = 1, and we 
still have the same point. We may take either direction 
along a line to be that of r increasing; if we reverse the 
direction the signs of all the are reversed for the same 
point. 


/ 
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For any point on (8), 

= Ukm («* + Im = ^ikm CcJm (25). 

This is a constant vector for all values of r, and may there- 
fore be considered as a property of the line. We denote it 
by If, Then we have six properties of the line given by 

l{j lt\ These are coordinates of the line. They are connected 
by two relations, 


= 1 (26), 

W = €ikm haJm = 0 (27). 

The 1/ have a geometrical interpretation. Thus if we con- 
sider the plane 

•*'2^3 ^3^2 ~ (28), 

this plane passes through the line. Also if — Cg > ^3 = > 

(28) is satisfied for all values of , and therefore if we take 

a line through a, parallel to the x^ axis, (28) represents the 

plane through ( 8 ) and this line. Two such planes determine 

the line, and therefore the and 1/ together determine the 
line. 

If we have two non-intersecting lines given by 


the line 


-f- rlf 

yt= Pi STTli 

rli 


passes through and is parallel to (29) 
(30) and (31) is, by (20), 


(29) . 

(30) , 

(31) 

The plane including 


^ikm ( 2 < - Pi) h'^m = 0 (32). 

This therefore represents a plane through (30) parallel to 
(29). The plane through (29) parallel to (30) is 

{Zi - a<) = 0 (33). 

The distance between these planes is the projection of the 
joining any two points on them upon a line perpen- 
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dicular to both. If a line perpendicular to both has direc- 
tion cosines , the shortest distance d between the lines 
is therefore given by 

d sin 6 = {pi — tti) rii sin 6 

=- ± iPi- at) 

~ i {^ikmPi^k'^m (34). 

^ikmPi'^m “ > ^ikm^i^k ~ (35)i 

and therefore 

d sin ^ = ± (- lifn/ - m,//), 

so that, apart from the ambiguity in sign, 

d sin 0 = Z,w/ + mil/ (36). 

Thus the shortest distance is directly expressible in terms 
of the coordinates of the two hnes. 

Now consider two intersecting hnes 

Xi = ai + Hi; yi = + smt (37). 

The area of the triangle formed by af^Xi, yt is ^rs sin 6. 
The projections of these points on the plane = 0 are 

(0, 02 , 03), (0, x.^j Xg), (0, 2/2* ^3) form a triangle whose 
area is 


1 1 

02 

<h 


1 

^2 

<h 

1 

^2 

^3 


0 



1 

y-z 

Vz 


0 

m^s 

m^s 


= irs {l^niy - Z3W2) 

= \rs7i^ sin 6 (38). 

Thus the projections of a triangle, and therefore of any 
plane area, on the coordinate planes are in the ratios of the 
direction cosines of the normal to the planes. A plane area 
can therefore be treated as a vector whose components are 
proportional to the direction cosines of the normal. 
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If we have a line given by 

Zi=^ Oi-h rli (39)^ 

and pfisa, point not on the line, let us suppose the point 
turned through an angle B about the Hne. The foot of the 
normal from to the line is given by 

r = /, - at) (40), 

and therefore the displacement from the foot of the normal 
to is equal to 

Pi - {Ui + I A (A - a,)) = (P, _ a.) _ 

= h {h (Pi — “i) — h (Pk — a*)} 

(41). 

The magnitude of this displacement, p, is given by 

P'^ = (Pi - - A (A - a,)Y 

= (Pi - «()“ - lih (Pi - a.) (Pk - (Zk) 

= i A (Pi - a,) - It (A - ak)}‘ (42), 

the J being needed because in the double summation each 
pair of values of the sufHxes would occur twice. 

The plane through and the hne is 

rti (zt - a,) = 0, 

subject to jiA = 0, 

^i (Pi - at) = 0, 

and is therefore 

(Zi - at) Ik (Pm - a„) = 0 (43), 

while the n, are proportional to £,i„4 (y3„ - a„,) and there- 
fore equal to ± (yS„ - a„)jp. 

If now we turn through an angle 0, it receives a dis- 
placement p {I — cos 6) along the normal to (39), and a 
displacement p sin 9 along the perpendicular to the plane 
(43). If it goes to we have therefore 

- A = - (1 - cos 6) Z* (p, - a,) - I, ip^ - a,)} 

± sin B €amh iPn, ~ (44). 
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If 0 is a small angle and we neglect 6^, the displacement is 
simply 

Vi- “ = =1= ^ (^m ~ ^m) 

and in particular 

yi - ^ l =± {?2^ (^3 - «3) - ^3^ (^2 - “2)} 

and from the additive form of tliis equation we see that the 
displacement is the sum of those given by separate small 
rotations 1^0 about axes through parallel to the co- 
ordinate axes. Conversely, displacements due to small 
rotations about axes through a point can be added vectori- 
ally as if all were apphed to the system in its original 
position, and give the same total displacement as if they 
were compounded into a single rotation about an axis by 
the vector rule. We still have, however, to establish a sign 
convention. We decide that 6 is to be taken positive if a 
turn about the axis of Xg is from Xi towards Xj. This would 
make 


Vi - A = - 0 (^2 - Oj) ; (^1 - Oi) (47) 

with = ^2 = 0, Z 3 = 1. Hence 


If we write 


Vi ^ikm^k^ i^tn ®m) 

1^6 ~ W/f , 

Yi ~~~ ^ikm^k “* ®m) 


(48) 


(49) 


For instance, if = 0, we have 


yi — ^1 — >'2 p2 ~ ~ 

ys ^3 ~ (^^)* 


For finite rotations we return to (44) or (41), keeping the 
positive sign in the second term. We may transfer the 
origin to a, to save writing. The coefficient of in 
is 

= (1 - cos ^) - sin 


(61)> 
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for & i; if A: = the coefficient is — (1 — cos 0) (1 — 

Thus 

^ik = (1 — cos 6) {IJj^ — S,-fc) - sin 0 (52), 

where c,^ is the antisymmetrical tensor corresponding to 
namely 

/ 0 ^ — ^ 2 \ 

(-?3 0 ^i) (53). 

\ h h 

Thus the displacement is represented in general by 
where is a tensor of the second order, expressed as the 
sum of symmetrical and antisymmetrical tensors. The anti- 
S5rmmetrical part is seen to be of the first order of magnitude 
in d and the symmetrical part of the second order. 


EXAMPLES 

1 * Given that the general quadric surface is 

S = + B^Xi + C = 0, 

prove that the locxis of the mid-points of parallel chords is a plane, 

and find the condition for this plane to be perpendicular to the 
chords. 

2. Find the condition that the line 

+ P = 0 


maj’ touch the quadric S. 


CHAPTER III 


PARTICLE DYNAMICS 

The essence of particle d3Tiamics is that the second de- 
rivatives of the coordinates of a particle with regard to the 
time are equal to functions of the position and velocity of 
the particle with reference to neighbouring particles. The 
relations therefore provide a set of differential equations 
to determine the coordinates. The equations for any 
particle can be put in the form 

mXi = 'LXi ( 1 ), 

where m is the mass of the particle, are the forces due 
to the other particles, and the summation is for all the 
other particles*. It is a matter of experiment that this 
form holds when the axes are a certain type of Cartesian 
axes, which we call non-rotating and unaccelerated, or in 
brief dynamical. 

If we take a different set of non-rotating axes at the 
same origin we have 

x/ = a^Xi ( 2 ), 

and since the axes are not rotating the are constants. 
Hence by differentiation 

( 3 ), 

x/ = aiiXi ( 4 ), 

and therefore velocity and acceleration are vectors. 

The force Xf on the particle due to some other particle 
is measured by the contribution to mXi due to the other 
particle; that is, the part of mXi that would disappear if 
the other particle was remov'cd to an infinite distance. It 

* Cf. Jeffreys, Scientific Inference, Chapter viii, for a fuller analysis. 
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is actually measured by the acceleration. If we have the 
three acceleration components of the particle 1 due to the 
particle 2, we can find the acceleration in the direction of 
by the formula (4), and the form (1) will become 

= SZ/ (5), 

provided that we define Z^' by the rule 

Z/ = a^jXi (6). 

The meaning of the force in any direction requires de- 
finition in any case: if we define it by its relation to the 
acceleration component in that direction, it follows auto- 
matically by (6) that force is a vector. It follows that the 
sum of any number of forces, obtained by adding the 
components separately, is also a vector. But the practical 
importance of the idea of force arises equally from the fact 
that in many cases the force components are known once 
for all from experience as functions of the coordinates and 
velocities. The form (6) then says, as a general principle, 
that the forces in any direction are additive. 

If we form the contracted or scalar product of mii, and 
£Z< by the vector we get 

mXiXi ~ TiXiXi (7), 

the left side of which is ^ {^mXi^). (The square of course 

implies the product of x^ by Xi and therefore the summation 
for i — 1, 2, 3.) We write 


T = {±1^ -h x^ -f x^) (8), 

and call T the kinetic energy of the particle. Then by in- 
tegration with regard to the time from to we get 




XiXidt 

to 



But in any short interval of time dt^ Xidt is the increase 
of w,, namely dx^. Hence the right side is equal to 

^ JC7 
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S j XidXi taken from the initial to the final position of the 

particle. We call the scalar Xfdxi the work done on the 
particle by the force in a small displacement dxi . Then 

XfdXi is the total work done by all the forces on the 

particle during the motion. We have therefore the scalar 
relation that in any motion of a particle 

Increase of kinetic energy = Work done on the particle 

( 10 ). 


In a system of particles we may add up the equations (10) 
for all the particles. We now take T, the kinetic energy of 
the whole system, as the sum of those of all the particles. 
Then we get 



the double summation implying summation for all pairs of 
particles; the first summation is for the particles producing 
the forces and the second for the particles acted on. It 
may happen that the Xt are all functions of the coordinates 
alone, and not of their velocities, and that provided the 
initial and final values of the Xi are the same the integral is 
the same however the Xi may vary in the interval. If so, 
the integral is the difference between the values of a certain 
function U for the initial and final positions of the system; 
we call the system conservative and U the ww‘k-function. 
Then (11) becomes 


or 


T ^ U = constant (13). 


This is the equation of conservation of energy. The quanti- 
ties on both sides are scalars. — U ia often denoted by V 
and called the potential energy. 
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If U exists, then in all possible small displacements of 
the system 

'Z'LXidXi^dU (14), 

where the 3n in number for n particles, are all in- 
dependent. The coefficient of dx^ for any one particle is 
the total force in the direction of Xt on that particle. 
Hence the force on a particle in the direction of x,- is 
dUjdXi, which is of course a vector. 

We can also take the scalar product of (1) by any set of 
small quantities Bx^ whatever that constitute a vector. 
Then 


mxtSXi = ^XiBxi ( 15 ) 

is a scalar relation; but as the are arbitrary we can 
equate their coefficients and regenerate the equation (1). 
If we now integrate with regard to t from Iq to we get 


rti 

mXiBxfdt 

JU 


1 “ 

JU 


'ZXiBXidt 


The left side, on integration by parts, gives 


(16). 



But we can consider the a:,- -f- Sr, as coordinates of a par- 
ticle in a motion differing slightly from the actual one ; that 
is, at a given value of t, the coordinates are r* -f- Sr^ instead 
of r< . Then 


d jj d , ^ . d 

dt ~~ dt 


dl 


= r< -h Sr, — Xi 

= (18), 

since Xi + Bxf is simply the varied velocity or rate of change 
of the varied coordinate Xf Sr^. Then 
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Then for every particle, to the first order, 

( 20 ). 

We can now add up these equations for all the particles of 

the system. If U exists we can express the result in the 
form 

8 (T + U) dt = 11, (21), 

<0 L J/o 

or, since the hmits are not varied, we can move the 8 out- 
side the integral. If the variation is such that the initial 
and final positions of the system are unaltered, 8a:< = 0 
when t = t^^QT and we have 

8 \^\t U)dt = 0 

JU 

to the first order in the variations of the path. This is 
Hamilton’s principle. 



{8 -f- 2Xi8x,}d^ = 


EXAMPLES 


1. If the Cartesian coordinates of every particle of a system are 
known functions of a set of generalized coordinates prove that 







dxf 



where the summation convention is also understood on the right 
side. Deduce Lagrange’s equations for the case where the q^. are 
all independent. 


2. The equations of motion of a particle are m£f = X,- — 
where k is constant. Prove that 

Hence show that for a system in periodic motion, or in one slowly 
changing its state, on an average over a long time, 2T = — 2x, X^. 



CHAPTER IV 


DYNAMICS OF RIGID BODIES 


A rigid body is one such that whenever it is displaced the 
distance between any two particles of it is unaltered. Since 
three particles By C are in a straight line if the sum of two 
of the distances AB, BC, GA is equal to the third, it follows 
that straight lines are unaltered by displacements of a rigid 
body. Since when AyByC are not in a straight line the angles 
of the triangle ABC are determinate functions of the three 
sides, it follows that all angles are imaltered by displace- 
ments of a rigid body. If three Hnes meeting at a point and 
fixed in the body are mutually perpendicular before dis- 
placement, they are stiU perpendicular after the displace- 
ment. 

The equations of dynamics in the form mXi = are 
true with respect to dynamical (that is, non-rotating and 
unaccelerated) axes. Let any particle Q of a body have 
coordinates Xt with reference to dynamical axes at 0. Then 
let the body be displaced in any way, and let the particle 
have the new coordinates a;/. We require the relation 
between a:/ and Xi . Consider a particle of the body, P say, 
whose coordinates before and after the displacement are 
and a/. Put 


= a, -f- ar/ = a/ + ?// (1). 

Then and y/ are the coordinates of Q with respect to 
axes at P parallel to the dynamical axes before and after 
the displacement. Also if we'tmagine the original axes at 
P to be specified by the particles on them, these particles 
in the new position still specify a set of rectangular axes, 
with respect to which the coordinates of Q are still If 
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the cosine of the angle between the y* axis in its new posi- 
tion and the y^ axis in its old one is , we have therefore 

Vi = (likyic ( 2 ), 

whence x/ = a/ + — a*) (3). 


The displacement of Q is 

a/ — (4). 

We can prove that there is a straight line of particles in the 
body such that the displacement of Q is the same as that 
of F. For such a particle we should have 


dikyk -yi= (<^ik - ^ik) yk = o (5), 

and these three homogeneous equations in are consistent 
provided 


that is, 



(6), 

( 7 ). 


®21 


®22 “ ^ 


a 
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Now the determinant || a,* || isxinity, and each element in 
it is equal to its first minor. If we expand (7) we get 

I I I I ~ {(^11^22 “ ®12®2l) + ~~ ®23^) 

4- (cEsaflii — UsiCEis)) + (cEii + ^22 d* ^ 33 ) “1 = 0 (8), 

since the terms all cancel. Hence (5) have an infinite 
number of solutions, all proportional, and the points there- 
fore lie on a straight line. 

If we take any two planes through this line, the angle 
between them is the same after displacement as before, 
and therefore all planes are rotated through the same angle. 
Thus any displacement of a rigid body is equivalent to a 
displacement of a particle of it combined with a rotation 
about a line through that particle. If the angle of rotation 
is dy we have, by comparing (6) with (52) of Chapter n, 

^ik — ^ik “* ^ik 


“ ihh - ^ik) (1 “ cos 0) - sin d (9), 
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where are the direction cosines of the axis of rotation. 
Thus the a**, are determined in terms of this axis and the 
angle of rotation. Conversely, if we take the symmetrical 
and antisymmetrical parts, 

{}ilk — S.-j.) (1 — cos 0) = ^ + 6 fc,) ( 10 ), 

J (cjtim “ ^ikm) Siu d — \ — 6jt,) (1 !)• 

But for a given w, with ikm in cyclic order, 

~ ^ikm = — 2 ( 12 ). 

Thus \ ( 6 ^* - 6 *^) (13), 

and 

sin^ e=i {( 6,2 - + (b^ - 632 )" + (^31 - 613 )^} (14). 

Now suppose that the displacement is small. Then a,* 
have nearly their values for zero displacement, that is, 
S^fc. The direction cosines of the and yi axes in their new 
positions are a,-fc and and thus, ii k ^ I, 

ai^aa^O (15), 

and, if A: = Z, = 1 (10). 

But in (15) for i ^ ky is nearly 1 and small ; for i = I, 
an is nearly 1 and afjg small; for i not equal to k or Z, both 
and an are small. Hence to the first order 

O'lk + Ufc( = 0 (17). 

If A; = Z, then for t # A; or Z both terms of ( 1 6) are small of the 
second order, and therefore for i = k or Z, o,;^ = 1 + a 
second order quantity. Thus to the first order 

^ik ~ — bile • ( 10 )» 

where 6 <* is an antisymmetrical tensor. 

The displacement of Q is Sa< -I- where Sa^ is 

the displacement of P, and is a first order small quantity. 
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Now let Xi be the total force acting on the particle at Q. 
In the displacement the work it does is the scalar 

hW = Xi (Sat + - Vi) 

~ “ ^ikyie) 

— XtSat — (Xi^ji + Xi^yt) — 16^ (Xtyjc — X^yt) 

(19). 

When i and k are interchanged, 6,* and Xtyj^ — Xj^yt are 
reversed in sign, while Xtyft + Xkyi\s unaltered. Hence 

Sir = Xthat + \ (ViXi^ — yjcXt) bn^ (20). 

Now write y.Z^ - y^Xt = Lt^ (21), 

so that Ltjf is an antisymmetrical tensor. The second term 
in (20) is the sum of nine terms, of which three are always 
zero and the others equal in pairs. If we replace 6^* and 
^ik by f-be associated vectors, we have 

SW = XMt + L^b^ (22), 

where the Sa, and b^ are all independent of one another 
and the same for all particles of the body. Hence if we add 
for all particles of the body 

Sir - (2Z,)Sa,+ (SLJ6„ (23), 

and the total work done in any given small displacement is 
determinate if we sum up the forces acting by the six 
expressions SZ^, Further, the contributions to 

SZf, 'LLfn from the internal reactions are zero. This follows 
at once if these reactions consist of equal and opposite 
forces between pairs of particles along the line joining 
them, and also has the justification that it leads to correct 
results. Then we may restrict 2Z,-, to the contribu- 
tions from the external forces. This is (TAlemberVs Prin- 
ciple. 

In the limiting case of continuous motion, we may con- 
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sider the displacements that take place in a short time 
interval Bty and put in the limit 

(24), 

and call Ui the velocity of P and a)ffc the angular velocity 
of the body. Then the velocity of Q is 


Xi — Ui — (^ikyje (25). 

Now consider the centroid (?, with coordinates Xi de- 
fined by 

(Sm) Xi ^ HrriXi (26). 


It is usually assumed without proof that G is fixed in the 
body, though this is not obvious. But suppose that the 
particle at O with coordinates Xi in the original position of 
the body goes to G' when the body is displaced, its new 
coordinates are x/, given by 


x/ — a/ "h (xjt — ttfc) (27), 

and the coordinates of the new centroid G" are given by 


(2m) Xi" = Unix/ (28). 

We have to show that G' and G" coincide. We have 


(2m) Xi — 2m {a/ -4- Oit (^a; “■ ®fc)} (29)> 

and therefore 

(2m) (xi x/) ^ 2m {a/ (Xf. — a*)} 

— (2m) {a/ + (r* — a*)} 

= 0 (30). 

Thus the particle originally at G is displaced to the new 
position of the centroid, and therefore the centroid is fixed 
in the body. 

Now return to the equations of motion of a particle of 
the body 

mXi = Xi (31). 
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By addition we form from these the equations 

'LmXi = HiXi (32), 

and by cross-multiplication followed by addition for all 
the particles of the body, 

Sm {XiZk - XkXi) = 2 (XiXk - Xj^Xt) (33). 


On the right sides of (32) and (33) the contributions from 
the internal reactions are zero, by d’Alembert’s principle. 
Also 



'LmXi = ^2 'LmXi = (2m) = Mx^ 

(34), 

where 

M = 2m 

(35), 

the total mass of the body. Also if 



hik = 2m {XiXk - x^Xi) 

(36), 

we can 

reduce (32) and (33) to 



Mxi = 2X( 

(37), 


i 

(38). 


These are the fundamental equations of rigid dynamics. 
The three independent hijg are expressible in terms of the 
three independent by (36), and we have therefore six 
differential equations for the and which determine 
them in terms of the initial conditions if the external forces 
are given. The motion of the body is therefore determinate. 

The principle of virtual work follows immediately. For 
if a body is initially at rest, J"i and are zero, and the 
condition that they may remain zero is that and 
shall vanish. But this imphes, by (23), that in any small 
displacement of the body the work done by the external 
forces is a small quantity of the second order in the dis- 
placements. Conversely, if there are six independent 
possible small displacements such that the external forces 
do no work, to the first order, in any of them, the co- 
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efficients of the independent and of (23) must all 
vanish. But these coefficients are the and and 
therefore the and A;* do not vary with the time. The 
vanishing of the work done by the external forces in six 
independent small displacements is therefore a necessary 
and sufficient condition for equilibrium. 

The equations (38) may be put in another form. If we 
consider any moving origin O', not necessarily fixed in the 
body, with coordinates a<, we can write for the coordinates 
of a particle with respect to O', 


Vi^Xi- at 


(39). 


If any vector associated with the particle, such as its 
velocity, momentum, or acceleration, or the force acting 
on it, has components we may form the antisymmetrical 
tensor yiUjg — y/fUt and call it the moment of the vector 
about O'. From (32) and (33) we can form the equations 

Sw (XiXn — Xj^Xt) — {at'Lmxjt — 


= 2 (XtXt, - x^Xt) - (a*2X* - a,.2X,) (40). 

where only the external forces make any contribution to 
the right side. But by (39) this is equivalent to 


2w (y^a-jfc - y^^Xi) = 2 (y.-X* - Vk^t) (41). 

Therefore the sum of the moments of the mass-accelera- 
tion products about any origin is equal to the sum of the 
moments of the external forces about that origin. 

If we denote the moment of momentum, or angular 
momentum, about O' by , we have 

Art' = 2m {yiXk — Vk^i) (^2), 

^ Art' = 2m (yiXk — yk^i) + - Vk^i) 

(43). 


and 
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The second term may be written 

2m {{Xf - Xk - (4 - df.) Xi} = - 2m (d.i* - dj^Xi) 

= -M(diX},^dkX,) (44). 

Thus (41) is equivalent to 

I h,,' + M (d,t; - = 2 {y,X, - y,X,) (45). 

In many important cases the second tensor on the left 
vanishes identically. This is clearly true if O' is fixed, when 
the time derivatives of the are zero ; when the centroid 
is fixed; and, if both are moving, if the vectors and 
are parallel, that is, if the velocity of the moving origin O' 
is parallel to that of the centroid. The most important case 
is where the moving origin is identical with the centroid, 
when the last condition holds automatically. These terms 
also disappear if the moving origin is an instantaneous 
centre of rotation always at the same distance from the 
centroid. They vanish for a sphere or circular cylinder 
rolling down an inclined plane, but not for a rolhng elliptic 
cylinder. 

If the moving origin is the centroid, 

h,f,' = 2m {yi (x* + y*) - Vk {^t + 2/.)} 

= 2m (y.y* - 2/*^,) (46), 

since 2my^ = 2my* = 0, 

by the definition of the centroid. Thus the angular mo- 
menta about the centroid are expressed completely in 
terms of positions and motions relative to the centroid, and 
the formulae for them have the same form as those for the 
total angular momenta with reference to a fixed origin. 

By (25) y,=^-w,^y^ (47), 

7/1 being here a dummy suffix. Substituting in (46) we have 

hik = - Sto - wt„ytyj 

= ~ ^Im<^k77, 


(48), 
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where is the symmetrical tensor defined by 

-Bim = (49), 

and depends only on the masses and positions of the par- 
ticles of the body. In the three-dimensional case, with i 
and k unequal, m must be equal in turn to i, k, and the 
other value. Then, for instance, 

— (-BaiWu + -B22W12 "I* 

“ (-Bn^a + Bj 2W22 + -B13CO23) 

~ (-^22 4 " -Bij) a»i2 -^ 23^31 (^ 9 ). 

If we replace the antisymmetrical tensors by the associated 
vectors this takes the form 


— (-Bu H" ■B22) W3 — B12OJ1 — ^ 23^2 


— (^11 -f ^22 + B^)c 02 — + B 23 OJ 2 + -BgjO^s) (51), 


and in general 



where is the symmetrical tensor given by 



- B..B 


kk^im 



( 52 ), 


= Si^ — 'Zmyiy„^ ( 53 ). 

It evidently corresponds to the system of moments and 
products of inertia given in ordinary treatises on dynamics ; 
-^11 > *^22 > -^33 Q’^e the ordinary moments of inertia Ay By C, 
but ^23 > -^31, ^12 are equal and opposite to the ordinary 
products of inertia F, O, H, 

It should be noticed that this reduction is characteristic 
of three dimensions; in a higher number of dimensions 
there is no analogous simplification of the form ( 48 ). 

The equations of motion of a rigid body then take the 
form 



( 54 ) , 

( 55 ) , 
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where Z-,*' is the moment of the external forces about the 
centroid. 

It may happen that one point of the body is fixed. In 
that case we may take the origin at that point. Then in 
addition to the known external forces there is a reaction 
at the origin which can if required be found from (37) ; but 
the reaction has no moment about the origin and the 


motion is given by (38). But in this case 

(56), 

and we find, by a process analogous to the last one, 

^tk = ^km<^im ~~ (57), 

where (58), 

In terms of the associated vector, 

(56), 

where <7,^ = S„„ - ILmXiX^ (CO). 


The correspond to the moments and products of inertia 
about the origin. They can be expressed in terms of those 
about the centroid ; for 

C,„ = (x* + y^Y) 8,„, - 2to (x, + j/,) (x„ + yj 

= 8 ,„ - 

+ (61), 

the terms linear in the y’s vanishing by the definition of the 
centroid. 

The relevant equations of motion then take the form 

(62), 

where tlie are the moments of the external forces about 
the origin, the reaction at the origin making no contri- 
bution. 
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The kinetio energy of the body is given by 


2T = 

= Sm {Xi - 

= + B^UiikOiim (63). 


In three dimensions is zero unless i is different from 

both k and m. If i = 1, A: = 2, m = 3, or if i = 1, A: = 3, 

m = 2, cojgcoig = — ioatog. If t = 1, A: = m = 2, 
and has coefficient B 22 * Thus o>^ enters with a coefficient 
^ 22 + ^ 11 » or But a> 2 tt >3 has a coefficient — 2^23 or 

+ ^ 32 . Thus in all 


2T = + Aj,m<^kOim ( 64 ). 

When a point of the body is kept fixed, 

2T = C,^w,oym ( 65 ). 


We notice that the linear and angular momenta take the 
forms 



dcDf 


{x„ constant), 



for origin at a point of the body held fixed. 

In one respect the forms (55) and (62) are inconvenient. 
They involve the tensors Af^n ^od 0^^? which depend on 
the Xf and i/< and therefore in general change as the body 
rotates. It is more convenient to use such axes that the 
relevant tensor in the particular problem is constant. To 
achieve this the axes must rotate, and then are no longer 
dynamical axes. Suppose then that we have a set of axes 
x/ rotating in any way, and that their direction cosines 
with respect to the dynamical axes x^ are ce,-^ . All the usual 
tensor relations hold for transformations from the 
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system to the Xi system and vice versa\ but the a,-^, while 
specified for each instant, are now functions of the time. 
If Ui is any vector, we have 

w/ = aiiUi\ Ui = aiju/ (67). 

Then - 37 ^ is another vector ; and 
at 



The component of this vector in the x{ direction is 


dut f duit 

ana,, + a,,u, 


(69). 


The first term of this is 


^ duj dui 

^'*~dt^'dt 


(70). 


Also 


is the Xi velocity of a point at unit distance along 
dt 

the x/ axis, that is, a point with coordinates (i = 1, 2, 3) 
with reference to the dynamical axes. This velocity is 
— dikdki^ where 6ik is the antis 3 Tnmetrical tensor expressing 
the rotation of the rigid frame consisting of the moving 
axes. Then 


anUj “tk^ki 


= - e,;u; 


(71), 


and the required component is 

dui 


dt 


- Bi/u; 


(72). 


If we use instead of di{ the associated vector, the three 
components become 

(ii/ — O 2 , Wg ^3^1+ Ui 6q , 

+ ( 73 ). 
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If for instance the components of a displacement along the 
a:/ axes are «/, the formula (73) gives the components of 
the velocity along these axes ; if are the components of 
velocity, (73) gives the components of acceleration; if 
are the components of angular momentum, (73) gives the 
components ol their rates of change with reference to 
dynamical axes, and these are equal to the components of 
the moments of the forces acting. 


EXAMPLE 


Prove that 





CH A PT ER V 


EQUIVALENCE OF SYSTEMS 

OF FORCES 

An external force Z; acting at a point a:,- of a rigid body 
produces dynamical effects summed up in the vector Z< 
and the antisymmetrical tensor ar^Zj^ — a^jtZ,-. A force has 
a line of action ; that is, if we take its resultant R given by 

= ( 1 ) 

we can define a direction If by 

Z, = Rh (2). 

By convention R is always taken positive. If a force Xt 
acts at the point a:^ + rlt , where r is arbitrary, we have 

{Zi H- rli) Zfc — {Xjg + rZfc) Z( = ^iXj^ ~~ x^Z^ (3). 

Thus the dynamical effects are the same if the force X{ 
acts at any point of a line through Xi with direction cosines 
proportional to Z^. This line is called the line of action of 
the force, and the force can be said to act along it. 

If a force has magnitude R and acts at ar^ in the direction 
Zj , we have 

Xi^Rli (4). 

Rik ~ R" (^)> 

or Ljn = Rim* (fi)) 

in the notation of the coordinates of a line. Thus the Z( 
and are the products of the resultant into the six 
coordinates of the line of action. 

The moment of the force Z< about a point a, is 

(Xi — a,) Zfc - (arjfc — a*) Z^ = i? (Z,V — a,/* + a^h) (7). 
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If we take the associated vector, its component in the 
direction m„ is 

li = R {fTlil/ + liVl/) (8), 

where m/ are the other coordinates of the line through 
in the direction . We notice that 

+ hm/ = d sin 6 (9), 

where 6 is the angle between the two lines and d the length 
of their common perpendicular. We may call Rd sin 6 the 
moment of the force about the line (m^, w/). 

By d’Alembert’s principle, the motion of a body is un- 
altered if to the forces acting on it we add two equal and 
opposite forces acting at the same point, or, by (3), along 
the same line of action. 

Now consider a pair of equal and opposite forces and 
— acting at points and They clearly make no 
contribution to XX f. Their contribution to is 

= {cii ~~ 6 ,) Xjf — (a/f — bfg) Xf ( 10 ). 

Since and are equally affected by any motion of the 

origin, the contribution of such a pair of forces applied to 

definite particles to both XX, and XL,,, is independent of 

t e position of the origin. Such a pair is called a couple, 

and its contribution to XL,„ is called the moment of the 
couple. 

If the vectors a, — b, and X, are both perpendicular to 
a line with direction cosines n , , the components of XL^ are 
t e components of a vector along this line ; while the forces 
act in the same plane perpendicular to the line. This plane 
IS called the plane of the couple, and the line an axis of the 
couple Evidently equal couples in parallel planes are 
equivalent. The magnitude of this vector is Rd, where 
« and d is the perpendicular distance of a, from a 
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line through 6,- parallel to and its components can 
accordingly be written Rdn„ . 

Any system of forces is equivalent to a force at an 
arbitrary point together with a couple. For with each 
Xi acting at we can associate a pair of forces ± at a<. 
Then the system is equivalent to EX,- at a,- together with 
a set of couples whose total moment is 

■^ik — 2 {{^i “ ®i) Xjg {xjg — a*) X,} (11)‘ 

But this is equivalent to a single couple ; for we have only 
to make proportional to the and Rd equal to their 
resultant. 

If are the moments about the origin, 

-Z/tt = 2 (XfXjt — XjgXi) (12), 

and if for brevity we replace 2X, by simply X^, 

~ ^ik ~ (^iXjf — OfcXd (13)> 

or = a, Xfc (14). 

Evidently X,’* is a scalar and independent of o,-. Also 

~ Ajt X„, (1^)* 

The first term is a scalar and independent of . The second 
is identically zero. Hence X,^ and XiM^ are scalar in- 
variants. 

The system is equivalent to a force X,- and a couple Mi 
at a^. These vectors are parallel if 

~ ^ikm^i^k ~ (16)i 

where p is a scalar length called the pitch. These give three 
linear relations between the three and the pitch, and we 
therefore expect a single infinity of solutions. But if we 
take the scalar product of (16) by X^ we have 

X^L^=pX^^ 


(17), 
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so that p is the same for all admissible values of . With 
this value of p, (16) represents three planes with a line in 
common. If we change to + aX ^ where a is a scalar, 
the left side of (16) is increased by 


“ ^ ( 1 ®)» 
and therefore if a, is one point on the Hne, all points on the 
line through parallel to satisfy the conditions. This 
line is the central axis of the system ; the system is equiva- 
lent to a force JH along the central axis and a couple G 
about it. This expresses the system as a torench. If we take 
the central axis as one of the coordinate axes, we have, 
since and X^Mf are scalars, 


(19), 

OR = XtMi - pXi^ (20), 

and therefore O, R^ and p are determined. 

The system can also be reduced to a couple parallel to 
a preassigned plane together with a force. For if S is the 
couple, and are the direction cosines of the normal to the 
plane, and if the force acts through a,, we have, for the 

moments about a^, 


= ( 21 ), 

and also ^ Sn^ (22). 

We have three equations to determine the and S. Again 
there are a single infinity of solutions. But if we take the 
scalar product by X„ , we get 

~ LfnXfn — 

= OR (23), 

so that 8 is determined provided n^Xm is not zero, that is, 
provided the resultant force is not parallel to the plane. 
Then the equations 


I'm ” 


( 24 ) 
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determine a line parallel to which is the line of action 
of the force. 

If li are the direction cosines of a hne through a,-, the 
moment of the system about this hne is 

~~ ( 25 ) 

= (26), 

where are the other coordinates of the line. If this 

moment vanishes the line is called a null line of the system. 
If bi is another point on it, 

Lm (K - flm) = €am<li^k {K - (27), 

which shows that hes in a definite plane through a,- . All 
null lines through a point therefore he in one plane. This 
plane is called the null plane of the point. 

Ah null lines in a plane pass through a point. For the 
system can in general be reduced to a couple in the plane 
and a force whose hne of action intersects the plane in one 
point. Then the system has no moment about any hne in 
the plane through this point, which is the null point of the 
plane. 

Any system is eqmvalent to two forces, one of which can 
be made to act along a given hne. For let the hnes of 
action pass through and 6< , and have direction cosines Z,- 
and and let the magnitudes of the forces be 8 and T. 
Then we have six equations, 

X, = Sli + Tm, (28), 

Lm = {aiSlt + biTm^) 

= Sl„' + Tm„' (29). 

The coordinates of the first hne being given, these are six 
hnear equations to determine the six coordinates of the 
second hne and 8 and T. But we have also 


= 1 ; niim/ — 0 


(30), 
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SO that we have in general just enough equations. Two such 
lines are conjugate lines. Clearly any line intersecting two 
conjugate lines is a nuU line. 

It can be shown easily that 

8 + Lilf) — LiXi (31), 

T2 = - 2SliXi + 8^ (32), 

whence (28) and (29) determine the coordinates of the 
second line explicitly. 


EXAMPLES 

1. A system of forces is reduced to a force at P together with 
a couple; P is chosen so that the couple is parallel to a given 
plane. Show that the locus of P is a straight line parallel to the 
central axis. 

2. Show that 

“m = 

is a point on the central axis. 

3. Two systems of forces are given by(X^, X,). Show 

that X^K^ + ^ is invariant. 


C H A PTER VI 


CONTINUOUS SYSTEMS 


In problems involving volume and surface integrals we find 
it convenient to denote elements of volume and surface by 
dr and dS respectively; both are always taken positive. 
They are of course scalars. The direction cosines of the 
normal to an element of surface are usually denoted by ; 
in most cases the normal is drawn outwards from the region 
under consideration. 

Green’s Lemma takes the form 



(1). 

and we have the corollary, if 


where F is a scalar, 

Ut — dV jdX( 

(2), 

11 


(3), 

where = 

02 02 02 02 

(4). 


and 0/0n denotes differentiation along the outward normal. 
Stokes’s Theorem takes the form 





= [j k^ikm dujdx^ dS 




The integral on the left is round a closed contour C. On 
the right is the direction cosine of the normal to any 
element dS of a surface S that fills up the contour; the 
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integral is over jS; and the sense of the normal is such that 
if the contour is described in the positive sense about any 
axis Xf , Z,- is taken positive when the normal is in the direc- 
tion of Xi increasing. 

The gravitational potential of a distribution of matter is 
given by 

V=fZy (6), 

where m is an element of mass, r is its distance from the 
point where V is to be found, and / is a constant equal to 
6*66 X 10“® when m and r are measured in grams and centi- 
metres and the time in seconds. When the mass is distributed 
continuously over a surface or through a volume, m must be 
replaced by adS or pdr respectively, where a and p are called 
the surface and volume densities. The work in displacing 
a mass m' through a small distance is m'dV* 

The electrostatic potential of a system of point charges 
is given by 

F=/S® (7). 

where e is a typical element of charge. If e is in electro- 
static units of charge and r in centimetres, /is + 1. The work 
in a small displacement of a charge e' in the field is — e'd V . 

The usual relations follow, that in free space in both 
cases 

= 0 (8), 

and in space occupied by matter of finite density 

V^F = - 477//> (9)* 

Also we have Gauss’s Theorem 



dS = — AmjYJm or — 


( 10 ), 


where the summation is over all the masses or charges 
within the closed surface S. In crossing a surface where 
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there is a finite surface density ct, dVjdn has a finite dis- 
continuity — 47r/o-. 

All these relations are scalar in form. For a gravitational 
field the force on a small particle of mass m' is 

m'Xi — m'dVjdXi (11), 

while for an electrostatic field the force on a small charge 
e' is 

e'Xi = - e'dVIdZi ( 12 ). 

The vectors Xi are called the intensities of gravitational 
and electric force respectively. The analogy in form be- 
tween (6) and (7) is constructed for mathematical con- 
venience; the difference in sign in (11) and (12) embodies 
the physical difference that whereas two positive masses 
attract, two positive charges repel. 

In a gravitating system we may construct a work- 
function 

(13), 

'VQ 

where the summation is over all pairs of particles m,, 
and Tpg is the distance between nij, and m,. The case where 
p = q IB excluded. Then the force on the pth particle is 
dWI{dXi)p . The potential at nip due to the other particles is 
Vp~ 1,'fmjrpgy the accent indicating that the case where 
q = p is excluded from the summation ; and the function 
"Znip Vp = 2Wy since each pair of particles is counted twice 
in this double summation. The function JSwp Vp therefore 
plays the part of a work-function. Similarly in electro- 
statics the function plays the part of the work- 

function with its sign reversed, that is, of a potential energy. 
These results may be generalized to the case of continuous 
distributions ; thus we can replace these functions by 

W = ^jjjpVdr + ijjcVdS (14), 
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where the first integral is through all regions where p is 
finite and the second over all surfaces where a is finite ; the 
first integral may therefore be taken through all space, 
excluding the surfaces where there is a surface density. It 
can be shown that with proper precautions about the de- 
finition of V the restriction p ^ q gives no trouble provided 
V is everywhere finite. 

Now if we consider the integral through aU space except 
thin laminae surrounding the surfaces where <j is finite. 



where the dn in the first integral is out from the region of 
integration and therefore towards the surface where cr is 
finite. On the two sides of such a surface the values of V 
differ by an indefinitely small amount, and for the two 


sides together 


dV^ 

dn 


dS = i^fadS, by Gauss’s Theorem. 


Hence the first integral in (15) is equal to — 



VdS 


taken over the surface and cancels the second term in (14). 
Hence 



through all space. In consequence of this form we may say 
that the gravitational work-function, or the electrostatic 
potential energy, is i^^/STr/per unit volume, where R is the 
resultant of the appropriate intensity vector. 

When the properties of the medium vary from place to 
place, V is no longer of the form S/e/r, and is no longer 
zero. But a potential still exists; if a small charge e' is 
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moved from a point P to a point Q, the work done is still 
the same whatever the route, and may be denoted by 
(Vp— Vq). The treatment is suggested by the fact that 
two similar condensers with the plates at the same po- 
tentials, but with air between the plates of one and another 
material between those of the other, have charges in a ratio 
K depending only on the media. dV/dn on the outside of 
the condenser being smaU compared with its value be- 
tween the plates, we infer that the charge per unit area for 
the same distribution of V is related to the discontinuity in 
KdV jdn, where K depends on the material. This suggests 
in turn that Gauss’s Theorem must be replaced by 



djv 

dn 


dS = — ^TrfJ^e 


(17), 


where the summation is for all charges inside S. Then 
applying this to the two sides of any surface we have 



and applying it to a region with a finite volume density we 
have 

- W III pdr = III A (Z |Z) dr (19), 
and as this must hold for all such regions, 



( 20 ). 


These equations are all homogeneous in F, />, < 7 ; hence 
the potential due to any set of charges is proportional to 
the charges if all are altered in the same ratio. Using this 
principle we can show by the usual method that the energy 
of a distribution is 




jlpVdr+ ^ jjaVdS 


( 21 ) 
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The first integral, applied to aU space except thin regions to 
cut off surfaces where there are surface densities, gives 


Stt/ 



V (k dr 


dXi \ dzj 


=-^ji 


KV^dS+ ^ 




dn ‘ Stt/ 

The surface integral cancels the integral \^^aVd8\ and 
hence 


so that the energy can be considered equal to KR^jSnf per 
unit volume. 

Magnetism may be treated similarly, starting from the 
assumption of volume and surface distributions of mag- 
netic pole strength, subject to the condition that the total 
pole strength in any solid is zero; or we may regard the 
ultimate magnetic unit as the doublet, which explains the 
need for the restriction involved in the former method of 
treatment. The potential at Zi due to a doublet of strength 
M at the origin with its axis in the direction A< is 


V = yMX,xj7^ (24), 

and if the doublet strength per unit volume in a solid is 1 
in the direction we can introduce the intensity of mag- 
netization at the vector = /A^, and say that the 
potential at is 

r = (25), 


where Ai corresponds to the point and dr = 
y is a constant. The magnetic force in free space is 
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We may write 




The potential is therefore equivalent to that due to a dis- 
tribution of magnetic poles liAt per unit area over the 

boundary and — dAijd^i per unit volume through the in- 
terior. 

Within a solid special treatment is needed. To define V 
or the force at , when Xi is within a solid, we must imagine 
a small cavity made about Xi , the intensity of magnetiza- 
tion everywhere remaining as before, and consider V and 

within it ; then the values of V and at x^ are defined to 
be the limits of those in the cavity when the dimensions of 
the cavity become indefinitely small. This process leads to 
little difficulty in gravitational and electrostatic problems, 
but in magnetism the limit of the force is found to depend 
on the shape and orientation of the cavity. The force in 
the cavity can be written 



(28), 


where V is given by (27); in the first integral the normal is 
inwards towards the cavity. The contributions to from 
the volume integral and the outer boundary are of the same 
form as for gravitation, and give no trouble. If the cavity 
is a cylinder in the direction of the intensity of magnetiza- 
tion, liAi is zero over the sides and equal to /, the resultant 
intensity of magnetization, on the ends. Such a surface 
density over the ends in the limit contributes nothing to V ; 
if the radius of the cylinder is small compared to its length 
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it also contributes nothing to but if the cylinder is 
of disc-like form it contributes 47ryA,/ to Hence if we 
take V for the complete body and define a< by 



(29), 


is the value taken by in a thin cylindrical cavity 
parallel to the intensity of magnetization. The force in a 
flat cylindrical cavity with its generators in this direction is 

a* = fl!( -t“ 4f7TyXiI = a,- -I- 4^yAi (30). 

Evidently and a,- are both vectors ; the former is called 
the magnetic force and the latter the magnetic induction. 
The theory of susceptibility and permeability may then be 
developed as usual. Also (27) shows that V is continuous 
across a boundary; but = dV/dn has a discontinuity 
— whence Zt-a,- is continuous across a boundary. 

The mutual potential energy of two doublets M and M* 
at and a;/, oriented in directions and A/, is 

W = \,'M' (31). 


where V is the potential at x/ due to the magnet at x ^ ; this 
gives 







yMM' 


(cos e — 3 cos B cos O') 


(33), 


where c is the angle between the axes of the magnets and 
0 and 0' are the angles made by the axes with the line 
joining the centres. 
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If the second magnet is turned through a small angle 8^ 
about a line with direction cosines 


by (49) of Chapter n. Hence, by (32), 




= yM3r ^ { a , 


3 cos 6 





so that the couple about a line parallel to the axis of x* is 



T 



Hydrostatics and Classical Hydrodynamics. The internal 
reaction in a fluid across an element of surface dS is a 
pressure pdS normal to that surface. If the density is 
the bodily force per unit mass Xi , and the velocity of the 
fluid at Xi is the acceleration of the fluid is found, by 
considering a small parallelepiped, to be given by 





If Ui is given in the Eulerian way as a function of the 
coordinates Xf and the time the operator djdt^ giving the 
rate of change of any element associated with a given 
particle of the fluid, is equivalent to 


a 

dt~dt'^^^ dx 



If we consider the circulation 12 around any closed 
circuit C in the fluid and moving with it, defined by 



( 3 ), 
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we have 


— * 


UidUi (4). 


The last integral is which always vanishes because 

when we move round the contour we come back to the 
same point, where the velocity has its original value. Also 
if Xi is the gradient of a single-valued potential, as when 
the bodily forces are due to gravity (the commonest case), 

I XidXi is the change of this potential round the contourand 

is zero. Again, if /> is a function of p only, as in an incom- 
pressible liquid or a gas at uniform temperature, the fii*st 
integral vanishes and 

If then O is ever zero around a circuit it remains so per- 
manently, This is true if the fluid is initially at rest and is 
set in motion by solids moving in it, and in various other 
cases of importance. But the vanishing of for all circuits 
is the condition for the existence of a velocity potential ^ 
such that 


d<t> 


( 6 ). 


In this case we can rewrite the equations of motion in 
the form 


dUi dui 


1 dp 

p dxi 




( 7 ), 


and mult^lying by dXi and adding we have 


0 9^, du 

it 0^/^' + 


-•‘dx,= --lldx, + X,dx, 

dXi p cXi 


smce 


dUi 


02 _ 0^ 
dXidXf^ dXi 


( 8 ), 


(9). 


JCT 
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This shows that for all contemporaneous variations 


^ = — f — H- 27 + constant 

J P 




where 


Z,= 


dJJ 

dXi 


( 10 ). 

( 11 ); 


is the square of the resultant velocity q. The constant of 
integration is not necessarily the same at all instants and 
therefore may be a function of the time. Hence we have 
the Bernoulli integral 




+ i? 


I^+U + F(t) 


(12). 


dt — j p 

The rate of change of mass within a given small paral- 
lelepiped dxidx 2 dx^ is equal and opposite to the rate of 
outflow ; hence g 

Z /fr =s — - — 


~ ai". 

and we have the equation of continuity 

^P _ d 
dt dxt 
dp dUi 

dt ^ dXi 


(13). 


or 


(14) . 

( 15 ) . 


Vectors with given Divergence and Curl. We sometimes 
have to find a vector w*- such that 

“bui 


dx, ~ ^ 
du^ 

dx.- dxu 


(1), 
( 2 ), 

where A is a given scalar and to,* a given antisymmetrical 
tensor. We want particular integrals of these equations. 

Evidently if 
where is any scalar, 


d4> 

- dx, 

(3), 

V^<j> = A 

(4), 

'-k 3“/ _ 0 

dxt 

(5), 
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and a solution of (4) is 




( 6 ), 


where ^ is to be evaluated at and A at f ^ ; 

dr = 

and r is the distance from to 


(7). 


= 


dF^_dj\ 

dXi. dx^ 


— €ikm 


dF, 

dx 


( 8 ), 


where is a vector such that 

dF^ 


dXf 


= 0 


we have 
Then also 


dUi 

dXi 


— ^i/cm 


= 0 


dUfi dUf dujg 

a^; 


dXidx^ 

dUu 


(9). 

(10). 


ikm 


d dF, 

“ 3^;. 0^p 

d^F, 

= — ^ikm ^Dfcs 

d^F 

= (8.»Sm»- 8.j.Sms) 0^^ 

. d^F, ^^F„ 

^‘•dXidx„ *’‘dxidx^ 

a dFt d^F„ 

“ 0x„ dXi dx,^ 

= - y^F„ 


by (9). 

Thus if 
we shall have 


F ^ 

477 


O^m 


dr 


( 11 ) 


( 12 ), 


du,f d\ii dUi ^ 


(13), 
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provided (9) is satisfied. But 


dF 


m 


dx 


m 


477 * 



CU 


m 


dx 


m 



- 1 dr 



(14), 

since is a function of alone and r a function of 
ii- -A^pplying Green’s Theorem to all space except a 
small sphere about we get 


dF 

dx. 


= — Lim 

477 - 


11 




dS ~ Lim 

477 - 


11 


f 1 dco 
r di, 


dr 

(15), 

since has the form of a gravitation potential and 
ml^^i that of a gravitational force. But the first integral 
vanishes in the limit when the sphere becomes very small, 
and the integrand in the second is zero provided the com- 
ponents Ut exist, for 


dco 


m 


_ d /dUf^ du,\ 
dx.^ \0a:< dxj 


dx 


m 


02% 


dxi dx 


= 0 (16). 


m 


If then we are given a scalar A and a vector such that 
its divergence is zero, a solution of (1) and (2) is 

dx^ 


u = + ^ 

* dXi dxk 


(17), 


where <f> and F^^ are given by (6) and (12) and ikm are 
in cyclic order. 

This analysis has two practical applications. In hydro- 
dynamics Ui is the velocity and is twice the vorticity, 
denoted by in Chapter ix. Here the divergence of the 
velocity and the vorticity may be given through all space, 
and the velocity (17) satisfies the conditions. If the actual 
velocity is v,-, we may put 

Vi= Uf-i- u/ 


dn/ 


cx 


= 0; 


du/ 
dXi dx,. 


= 0 


and then 


(18). 

(19). 
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Thus Ui is the gradient of a scalar <}>' satisfying Laplace’s 
equation. There is no such scalar that makes the velocity 
finite everywhere, including at an infinite distance, except 
a linear function of the coordinates. Hence Ui is the same 
everywhere. 

If there are solid boundaries or free surfaces at a finite 
distance may not be constant. If the region where 
vorticity is present does not extend to a boundary, there is 
no contribution to (14) from points outside this region; we 
therefore take (14) through a region large enough to contain 
the whole of the vorticity. Then the surface integral in (15) 
must also be taken over the boundary of this region, but 
still vanishes, and (17) is still a solution. But (17) may not 
satisfy the boundary conditions, and then we must add to 
Ui an irrotational solution chosen to make the whole velocity 
satisfy them. 

In electromagnetism Ui may be the magnetic intensity 
and the electric current across unit surface in a plane 
of constant. 

In many cases Ui has no curl outside a limited region of 
very small cross-section. This is often true in the motion of 
a real fluid, when the region may be called a vortex fila- 
ment, and in magnetism, when the region is a wire carrying 
an electric current. The former statement may be expressed 
also by saying that the motion is irrotational outside the 
vortex filament ; the latter says that magnetic forces due to 
electric currents have a potential. In either case the integral 

= [ UidXi (20), 

Jo 

taken around a closed circuit, is zero if the circuit can be 
filled up by a surface not cutting the filament or the wire, 
and has a constant value for all circuits that cannot be so 
filled up. The two conditions are mutually exclusive, and 
therefore the critical region itself must be a closed circuit. 
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In each case A vanishes. Then ^ = 0, and the component 
F„ is given by 

taken through the critical region. But if we consider the 
contribution from an element between two planes sepa- 
rated by dfm, and call the element of surface in a plane 
parallel to these dS^ we have 

dr^di^dS (22), 

fJoj^dS = = n (23), 

where the line integral is taken around the boundary of 
the filament. Hence 

- B It' 

the integral being taken around the length of the filament. 
Also 

47r (e^jt Jr dz^J r ) 

- a 1 ai © 


_ Q f 

47r f 


^ikm 


(it - Im 


ds 



where ds is an element of length of the filament and a 
direction cosine. 

In hydrodynamics H is the circulation around the fila- 
ment. In electromagnetism the unit current is such that 
if it flows in a circle of radius 1 cm. it produces magnetic 
force 277 at the centre. If we take the circle to be in the plane 
of constant, with its centre at the origin, we have 

ds = do \ 7 = 1 ; aTj = = a :3 = 0; 

= cos 0, ^2 = sin Oy ^3 = 0; = — sin 6, = cos 
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and 

j ^ = 0 for i = 1 or 2 

and = 27t for i = 3. 

Also for unit current , the magnetic force, = 277 . 

Hence, in this case, by (25), 

H = 477, 

and in general, if the current is /, Q = 477 /, and 

= I \ ds (26) 

= 1 1 dU (27). 

This may be transformed by Stokes’s Theorem into an 
integral over a surface with the wire as its boundary ; thus 

Ui I dS 

= - 7 (-) dS 




(28). 


But — - ■ (-\ is the magnetic force at Xt due to 

\^ / 

a doublet of unit strength at with its axis in the direction 
Im' Hence the force at any point is equivalent to that due 
to a distribution of doublets over the closing surface, with 


intensity / per unit area and directed normally to the sur- 
face. Such a distribution constitutes a magnetic shell of 


strength I. 
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It is customary to assert as a fundamental postulate that 
the magnetic force due to an electric current is eqmvalent 
to that of a magnetic shell, but I think this course un- 
desirable. The magnetic shell does not exist in nature and 
direct experimental test is therefore impossible. Further, 
though it gives the same force, it does not give the same 
potential ; the magnetic potential due to a magnetic shell 
is a single-valued function with a finite discontinuity at the 
shell, while that due to an electric current is a cyclic 
function with no discontinuity except at the wire itself. 

Mutual Energy of Electric Circuit and Magnetic Field. 
A magnetic pole m at a:, is under a force mUi ; but 


mUi — — ml 


dX( jj 


I 


m 



(29) 


and — l^ is the magnetic force normal to the 



surface at due to the pole. Hence - \\ljn^ 

is the total flux across the closing surface, or through the 
circuit, of the magnetic force due to the pole. Also if F is 
the magnetic potential due to the current, 


dV 

U4= — 


dXi 

and therefore the correct forces are given by taking 


(30), 


V = -l\^Nd8 (31), 

where N is the normal magnetic force due to a unit mag- 
netic pole. This can now be generalized; we say that in 

general the mutual potential W of a current and a set of 
magnetic poles is 


W = - / 


NdS 


(32), 
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where N is now the normal magnetic force due to the whole 
of the rest of the field. This form may be further extended 
to express the mutual energy of two electric currents. 
From this equation we can infer, as is done in the standard 

works of Jeans and Livens, that if jjNdS varies with the 

time, the variation generates an induced e.m.f. — ~ IJncIjS 
in the circuit. 

If we have two circuits carrying currents 7 and J, their 
mutual influence is expressed by the statement that their 
mutual energy is of the form (32), where N is taken to be 
the force near f ^ due to the current J ; points on the latter 
circuit can be taken to be given by . Then 

W = IJM (33), 


where 


M 




dS 


(34), 


and Ui is the magnetic force at due to unit current in the 
second circuit. But if Fi is the vector potential due to such 
a current, 

dF., 


M 




BF. 


= \F^d^ 


m 


(35), 


(36) 


(37), 


taken around the circuit. But 


F 


m 




dx 


m 


(38), 


taken around the second circuit; and therefore 

dx^dS„ 


M 


- 


(30), 


taken around both circuits. This gives the required form 
for the coefficient of mutual induction of two circuits. 



CHAPTER VII 

ISOTROPIC TENSORS 

A tensor is called isotropic if its components retain the 
same values however the axes are rotated. We have 
already obtained three examples, namely 8,i, and 

^iks ^mv8 • 

There are no isotropic tensors of the first order. For if 

M. was such a tensor, let us give the axes a small rotation 

expressed by the antisymmetrical tensor c,*. Then in the 
new system 

< == (S,/ - Cii) Ui = Ui~ CiiUf (1), 

and this can be equal to Uf only if 

CiiUf = 0 (2) 

for all admissible values of the . Thus 

CiiWj + C12U2 + C23W3 = 0 

C21W1 + + C23W3 = 0 ( 3 ). 

C 3 iWi + C32W2 + C33W3 = 0 

ButCjj = C22 = C33 = 0, while C12, 0^3 , Cgj are independent and 
equal and opposite to the components obtained by inter- 
changing suffixes. Hence ( 3 ) can be satisfied only if 

= ^2 = W3 = 0 ( 4 ), 

and therefore there is no isotropic tensor of the first order 
other than zero. 

If w,jt is an isotropic tensor of the second order, 

— c„) — Cjti) Uji 

= Uijg — CijhjfiUfi — 

~ "^ik (1) 

to the first order, for all values of i and k. Hence 


( 2 ). 
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If % and Jc are unequal, take i = 1, A: = 2. Since — 0 

we have 

<02^^ + <a3^2 “1“ ^"*^3 ~ ^ (3)> 

and therefore 

^ = 0 * ^11 = ^22 ( 4 )* 

By symmetry Ui^ is therefore 0 if i k, while = ^422 = "“33 • 

If i and k are both 1, we have 

H" ^3^31 H" ^'*4x2 “f" ^3 '*^3 ~ ^ (^)i 


which is satisfied since every term vanishes. Hence the 
only isotropic tensor of order 2 is a scalar multiple of S.-^. 
If Uifgfn is an isotropic tensor of the third order, 

ifcm ~ (^A:* ^fcl) (^mn ^mn) ^jln (1)> 

and therefore, for all values of i, k, m, 

^ ( 2 ). 

Take i = k = 1. Then 

"1“ ^13'*%lfn ^12^12m "1" 4'i3^13Tn 

"f“ Cfnl'*4iii + C,„2Wj22 ^ (3). 

Now put m = 2 so that 0^2 = 0- Then 

W212 + ^22 = ^111 

W312 + '^132 “ '*^12 ' (^)* 

'*4ii 3 fi 

From the last equation, and by symmetry, = 0 if two 
of i. A:, m are equal and the third unequal. Then by the 
first, is also zero if all of i. A:, m are equal; and the 
second shows that 

"^ikm — '^kim • 

If in (3) we put m = 1, every term vanishes, so that (3) 
holds. 

Now in (2) if i. A:, m are all different, is zero unless 
j = i, and then = 0. Hence (2) holds. It follows that the 
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only isotropic tensors of order 3 are scalar multiples 

of fjfcm* 

If '^ikmv is an isotropic tensor of order 4, we have, 
similarly, 

^mi'^ikip + ^pj'^ikmi ~ ^ ( 1 ). 

There are only three possible values for k, rrij p, and 
therefore at least two of them must be equal. We may con- 
sider separately the cases where (a) two are equal and the 
other two unequal, (6) three equal, (c) two equal and the 
other two equal, (d) all four equal. 

In case (a), take i = k = 1, m = 2, p = 3. Then 

^12^2123 “1* ^13^3123 "b ^"^223 "b ^13^1323 

-b ^21^1113 + <^23 ^33 "b C3iWji2i + C32 14^22 = ^ {2)» 

Hence, by the antisymmetrical property of 

^2123 + ^223 “ ^113 = ® 

^123 + ^1323 “ '*^121 = 

^33 “ ^122 = 0 

Other instances of case (a) can be obtained by inter- 
changing suffixes that are not already equal, and by turning 
the axes so as to bring 3 into the position of 1, 1 into that 
of 2, and 2 into that of 3. Thus (4) gives 

^1133 = ^1122 ~ ^2233 = ^2211 “ ^3322 “ ^3311 (^)* 

And also 

^1313 = ^1212 ” ^2323 ~ ^2121 = ^3232 “ ^3131 (6)» 

^3113 == ^2112 = ^3223 = ^X221 ~ ^2332 = ^1331 (^)' 

In case (6), take i = k = m = l,p = 2. 

^12^2112 “b CJ3U3JJ2 + C12W1212 + ^13^^^1312 

"b ^12^1122 "b Cj3Wjj32 -f- + ^23^^113 ~ (^)* 

The last term shows that 


(3) , 

(4) . 



^hll3 = 0 
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and therefore, by interchange of suffixes, all components of 
class (6) are zero. Also, from the coefficient of Cjg, 

^112 *!■ ^1312 "^1132 ~ ^ (f^)* 

But in (3) the last term vanishes and we infer 



^2123 + “^223 “ ^ 

(11). 

and therefore 


'^1312 "h '*^ 3112 “ ^ 

(12), 

whence, by (10), 


^132 ~ ^ 

(13). 


Thus all components of class (a) are also zero. 
The coefficient of Cjg in (8) gives 


^im ~ ^2112 "I" '*^1212 "J" '^1122 (^'^)» 

so that the components of class (d) are expressible in terms 
of the three types of class (c). 

No further information is got by transforming com- 
ponents of classes (c) and (rf). Thus if A: = l,m = p = 2, 
replacing i ox k hy j will give a zero component unless^* is 
equal to 1 ; and then the factor or is zero, and the 
relation holds automatically. Similar considerations apply 
if all of 2 , k, TO, p are equal. 

We may denote the components of type (5) by A, those 
of type (6) hy + V, and those of type (7) hy p — v. Then 
(14) gives 

Wjll] = '^2222 “ '^3333 — A -f* 2/X (15). 

There appear therefore to be three independent isotropic 
tensors of order 4, obtained by taking each of A, p, v in 
turn equal to 1 and the others to zero. 

In the A tensor, = 1 if t = and m = p, and in all 

other cases is zero. It is therefore equivalent to 
which is obviously a tensor of order 4, being the product of 
two tensors of order 2. 

In the p tensor, = 1 ^ i = to, A; = p, or if i = p, 
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Ic = nij and k. H also i = the component is 2. Other 
components are zero. This can be written 

and is obviously a tensor of order 4. 

In the V tensor, — 1 if i ~ mj k = pt and = — 1 if 
i = p, k — m, and in all other cases is zero. If also i — k, 
'^ikmp is zero. In this case, therefore, 

This can also be written 

'^ikmv ~ ^iik^miv (1S)» 

for if i = 1, = 3, = 0 unless J = 2 and then = 1. But 

then =1 ifwi=l, ^>=3, — 1 ifw=3, p=l, and 
otherwise — 0. Thus 


^313 “ ^331 ^ 

with corresponding values for the other components. 
Evidently (17) and (18) represent a tensor of order 4. It 
has already appeared in Chapters i and vi. 

The general isotropic tensor of order 4 is therefore 

SipSjttn) + »' (Sim^fcp” (2®)» 

where A, /x, v are scalars. 


EXAMPLE 

(S,mSfrp + S(pSjtm) "^ik = ^mp + ^pm* 

W,fc =* ‘^mp ~ ^pm- 


Prove that 
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ELASTICITY 


In an elastic solid, as in a fluid, the distance between any 
two particles of the body usually varies with the time. 
The body, however, has an equilibrium configuration that 
could persist if the external forces were zero or constant. 
We may take this as a standard of reference. If a 
particle actually at would be at x,- — in the standard 
configuration, we call Ui the displacement at x, ; in practice 
the squares of the Ui can usually be neglected. Evidently 

is a vector. 

At a point x< + where yi is small, the displacement is 
Ui-\- Vi, where 



( 1 ) 

( 2 ), 


where c,*. and are the symmetrical and antisymmetrical 
tensors 


^ik 


, (dut , 0mA , , 


If firt is zero at x^, the displacement has the same form 
as that due to a general displacement Ui together with a 
rotation expressed by the tensor Also, if is every- 
where zero, 




( 4 ). 
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Hence the are constants and the rotation is the same 
everywhere. The vanishing of the c.-^ is therefore the neces- 
sary and sufficient condition for a pure rotation. 

Now consider the part of due to the e,*. If we consider 
the quadric surface 

^ikViVk = ( 5 ), 

where r is a constant chosen so as to make the surface pass 
through the normal at y,- has direction cosines pro- 
portional to Cii-yff and therefore to The displacement 
due to c,t is therefore parallel to the normal at y< to this 
quadric. This coincides with the direction of the radius 
vector to y, if y, is on a principal axis. There are therefore 
three directions such that the relative displacement due 
to the is in the direction of y^, and these directions are 
mutually perpendicular. If we take new axes ^2» ^3 
along them, the quadric reduces to 

( 6 ). 

and all terms with j ^ I are zero. The displacement in 
the direction is now , so that aU distances in this 
direction are increased in the ratio 1 to 1 -f tyi'. The dis- 
placement due to the e,fc is then the resultant of three 
homogeneous strains parallel to three orthogonal axes. 

We see therefore that the displacement in any small 
neighbourhood can be represented as the combination of 
a rotation with three extensions at right angles. The latter 
express the changes of size and shape of an element of the 
solid. For this reason 6,^ is called the strain tensor and 
the rotation tensor. Evidently e.t has six independent 
components. For we can have 

= eyj, v^ = v^ = 0 (7), 

making all the c,*. zero except Cn = e; similarly C22 and ^33 
may exist independently of the others. Also, if 

= 0, ^2 = ^3 = 


( 8 ), 
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all the Cif. vanish except e^2 ~ Similarly €31 and 612 
can be assigned independently. 

In an elastic solid the internal force across an element of 
surface is in general inclined to the surface. If the area of 
the element is dSy the force across it must be specified by 
three components of the form PnidS parallel to the axes; 
where n may be regarded as indicating the normal to the 
surface. If we consider a small parallelepiped with edges 
dx ^, , dx ^ , dx^ y centred at {xi , x ^ , 0:3), the force across the face 
of area dz^dx^ centred at z^ + \dz^ is (p„, Pia, Pu) dz^dx^y 
where the are evaluated at aji + \dz^. The force on the 
opposite face is — (Pn, pi2, Piz) evaluated at 


Xj — Jdxi, and the total is dz^dx^dx^. In general the 
force in the z^ direction due to the stress across the faces 
of Xjfc constant is dz^dx^dx^y and we take account of 

As 

aU faces by using the summation convention. If the 
acceleration of the element has components and the 
density is p, the mass is pdx^dx^dx^-, while if the bodily 
force acting has components Xi per unit mass, the equa- 
tions of motion are 



The system of quantities p.^ constitutes a symmetrical 
tensor of the second order. To prove this, consider first a 
plane whose normal has direction cosines , intersecting 
lines through x< parallel to the coordinate axes at short 
distances from Xi\ thus a small tetrahedron is formed, 
whose sides are of order ly say. Let the area of the sloping 
face be dS\ then those of the others are a^dS. Consider 
now the rate of change of momentum of the matter within 
this tetrahedron. Evidently the contributions from the ac- 
celeration and Xi are of the order of the volume, that is, 
of Z®. The force across dS has magnitude p^dS. That across 
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the face of Xj, constant is — p^i times the area of the 

face, that is, — PkiO^ki^^’ But dS is of order Hence we 
have 

(Pii - aksPki) 0 (l^) = 0 (/3) (10), 

and hence if I is indefinitely small we have for the stress 
across a plane normal to at Xi , 

Vii^o,k}Pki (11). 

Now consider three perpendicular directions with direction 
cosines {j = 1, 2, 3). The force per unit area across a 
plane perpendicular to one of these axes, in the direction 
of Xi , is given by p^i . Resolving this along the direction of 
one of the new axes Xi we have, therefore, 

Pn — ^ki^ixPki ~ ^ii^klPik (12) 

by interchanging i and k\ which is precisely the law of 
transformation of a second order tensor. 

Consider again a small parallelepiped centred at with 
edges parallel to the coordinate axes, and form the equa- 
tion of moments about its centre. The contributions from 
fk and Xjf are of order at most, where the edges have 
lengths of order 1. The moment about an axis parallel to x^ 
of the stress p^^ in the face x^ — constant and parallel to Xi 
is the product of p^i into the area of the face and the dis- 
tance of the face from the centre; that is, to order i®, 
\P 2 \dx^(lx 2 dx^ . The opposite face makes an equal contribu- 
tion. The stress pjg in the face x^ = constant and parallel 
to tends to turn in the opposite direction. We have, 

I n 

{Pz\ - Pn) dx^dx^dxs = O {I*) (13), 

and therefore when we make I indefinitely small we must 
have 

Pzi ~ Pit (1^)» 

and in general 

Pik =" Pki (15), 

so that the tensor p^^ is symmetrical. 
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Now consider the energy interchange between the small 
parallelepiped and its surroundings. The stresses across the 
face of area dx^dx^ centred at + \dx^ are doing work on 
the element at a rate (ptiu*) dx^dx^t and the contribution 
from the two opposite faces is 



dx^dx^dx^ 


Thus in all the stresses are doing work at a rate 



(Pki<u) dr. 


The external forces are doing work at a rate The 

kinetic energy of the element is and is increasing 

at a rate pUiffdr. (We consider the actual specimen of 
matter occupying the element of volume dr at time t ; thus 
its mass is pdr, and if we keep to the same piece of matter 
at time t dt the mass is unaltered. If we considered the 
variation of energy within a given element of volume we 
should have to allow for the variation of p and the fact that 
the matter moving out of the element is taking its energy 
with it.) The rate of performance of work on the element 
therefore exceeds the rate of increase of kinetic energy by 


f a 


[a^ + p^i^i — p^ifi 


• dr 


(_a 

dx 





by the equations of motion, 


= ( 16 ). 

This work is stored up as internal energy of the element of 
volume. Evidently from its form it is a scalar. 

In any elastic solid the internal energy is a definite 
function of the state of the solid. In any change of state 
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the increase of internal energy therefore depends only on 
the initial and final states and not on the method of passage 
from one to the other. Now we have seen that six of the 
e,fc are independent, and if the element acqxiires displace- 
ments hui in time ht the corresponding increase of internal 


energy is 


,-8 dr. There is an apparent asymmetry 


\dxj 


according as i and k are equal or unequal. Thus p^dr has 
coefficient S (dv^jdx-^ — Scu, but i = 1, k — 2 contributes 
^21 S (dujdx^, and i = 2, k = 1 contributes (du^jdx-^), 
the two together giving 2j^j2 8€i2. But this is the same as 
Pn S€i 2 4- P21 ^^21 > whole contribution from the 

changes of strain is pijgZeijgdr. Also during the process an 
amount of heat hQdr may be absorbed. If then Edr is the 
internal energy of the element, 


hE — Pik^^ik 4" (l"^)* 


Since ^ is a definite function of the state of the system, 
and six of the 86,* are independent and determine the 
other three, 8^ depends on the changes of the temperature 
and of the six independent €,fc and has a definite value in 
whatever order these changes come about. But 



If the absolute temperature is 0, and a certain amount of 
heat 8Q is absorbed without change of any linear dimen- 
sion, the rise of temperature is related to 8Q by the rule 


hQ = pche (19), 

where c is the specific heat at constant strain. If there is 
also a change of strain, since 8Q and 80 are scalars, we must 
have 

= Qik^^ik 4- pcB6 (20), 

where the constitute a tensor of the second order. But BE 
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and SO/5 are perfect difiEerentials. Hence if we replace the 
six independent by we can write 



BE = {p^ + q^) Se,. + pc 85 

(21), 



(22), 

and 


g 

ai; + 5*) = 

9 9 9 

ai; : se ^P' + = 07, 

(23). 



It follows at once that if 5 is kept constant 
2g,8c, are perfect differentials. Also 


de ~ ^ ai; 







(25), 


and therefore 

g. = - e (26). 

If 85 = 0 we can write 


where 


£p,8e,= 8IF (27), 

2 IF = Co i- 2c,e,. + c„c,c, + 0(c®) (28). 


The Co, c,., c„ may be functions of 5. Then if we retain 
only terms in TF up to order e=*, 


Pr= Cr-^ Crse, (29). 

The Cf represent the stresses that would remain if the 
strains e, were removed without change of temperature. 
In most practical cases the original state is one of uniform 
temperature and no stress, so that c,. = 0. If there is a rise 
of temperature 6' under no stress, an element will acquire 
displacements 

(30), 


Vi “ aikO'yk 
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where a,* is a second order tensor expressing the thermal 
expansion. Thus 

^ik — \ («ifc "f Cfci) (31)i 

where is a symmetrical tensor; and 

Pik — Cifc + Cikmppmp^' ( 32 ), 

where „,p is a fourth order tensor. 

But by h3^pothesis this deformation takes place under no 
stress and therefore put = 0. This determines c,fc, and our 
formula for the stress is 

pik — ^ik,mp i^mp ~ Pmp^ ) (33). 

The coefficient c,fc„,p is the coefficient of in IF. 

Since there are six independent e,. there are twenty-one 
possible terms in a quadratic form c„e,.ea, and therefore 
there are twenty-one coefficients c,* • They clearly form 

a tensor of order 4; such a tensor in general would have 
eighty-one components, but this satisfies the symmetry 
relations that it is unaltered if we interchange i and k, or 
m and p, or i and k together with m and p together. 

From (2G) and (33), 

0 0 

“ ^^mp 0^ (^i*,mp) "1" ^ 0^ (^it, mpPmp ^ ) 

0 

~ ^^ik,mpPmp (34) 

if 6' is small. The second term does not involve the 6^.; 
the first is small of the first order in the e,*. 

IMany solids are isotropic; that is, they have the same 
properties in all directions. This applies to vitreous (glassy) 
sohds and to mixtures of crystals oriented at random. In 
that case a uniform rise of temperature in an element 
gives an equal expansion in all directions and 

V( == aO'yt 


(35) 
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simply; then a is the coefficient of linear expansion and 

= aSrt (36). 

The second order terms in constitute a scalar; and 
we have 


- C 


(37), 


a tensor of order 4. If it is isotropic it must be of the 
form (20) of Chapter vrr. Then the linear terms in give 

'Pik ~ ( 38 ) 

= AS,fcS„j,emp -f fj. + S,p 8j^) 

+ >' (3,m 3jtp — 8,p €„p 

— AS, + /i + SipCfcp) + V 3,pCjtp) 

= A8rtC„„ + /i (e^ + Cfc() + V (Cik - e*,) 

= AStte„„ + 2yxCrt (39), 

the last term vanishing since e,^ is symmetrical. 

Then 


^^2 = Pik^tk (40) 

= A A* + 2/i {e„2 + £^^2 + £332 + 26232 + 2e3i2 + 2ei22) (41), 

A — = dUifdXi (42) , 

The scalars A and fi represent properties of the material. 
Both are positive. 

VVe can also write 

21^2 « (A + 2/a) a* + 4/a (£23* + ^31^ + «12* 


^^11 ^11^22) (^ 3 ), 

This appears to differ from the form in Love’s Elasticity^ 
1906, p. 100, but the present £,* differ from Love’s strain 
components. My is the same as his, namely dujdx -^ ; but 
icy ^23 is only half his, so that his assemblage of strain 
components is not a tensor. 
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If all the Cii were equal to one another and therefore to 
|A, we should have 

21^2= (A+ f^)A2 (44), 

Pik = (A + A (i = k ) ; =0 k) (45). 

In general we write 

A + §/x - (46), 

and call k the bulk-modulus. 


2 irj=(A 

+ 

Ip) A“ 






+ 

V (eii^ 

+ 


633^ + 2 ^ 23 ^ 

+ 263,’“ + 2e./) 



Ip 

+ 


633^ + 2 c 22®33 "i" 2633 4 " 2622622) 

= (A 

+ 







+ 

3 ^ (^11^ 

+ 


633^ ~ 

— ^611 “ ®11^) 



4/a 

+ 


ex2^) 


= (A 

+ 







+ 

{(^ 

— 

^33)^ "f" 

(633 — 621)® 

1 

+ 


+ 6623^ + + 6Cj2^} (47). 


The coefficient of /i vanishes if and only if the strain is 
a symmetrical expansion, and may therefore be called the 

distortional strain energy. 

If we allow for variations in temperature, 


Vik = A (A — 3ce0^) + 2/a (e,jt — ad 5,*) 

= {AA - (3A + 2/a) ad') 8,^ + (48) 

and qijc = — 0 ?q (49)- 


Every term in p.t is of the first order in the displacements; 
but dd' jdd — 1 and therefore gives rise to a constant term. 
This term is 

g-.-jt = 0.3A:£eS,jfc 


(50). 
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If 8Q = 0, so that no heat is lost or gained by conduction, 

pchd = — Qik^^ik ~ — 3/; ad 

= - 3^fcad8A (51), 

and therefore, if the strain takes place adiabatically. 


and 


d' = _ A 
pc 

/v 9A;^a*d\ ,v . «. 

Pik = H — j S<fcA + 


(52), 


(53). 


Thus in an adiabatic disturbance the constant A is in- 
creased above its value for a standard disturbance to A', 


where 


A' - A-}- 


pc 


(54), 


while pL is unaltered. The bulk-modulus k is therefore also 
increased to k\ where 

k' = k + 5^^ (55). 

pc 

In a simple thermal expansion at zero stress the absorp- 
tion of heat SQ is equal to pc^Sd, where Cp is called the 
specific heat at zero stress, and is the specific heat mea- 
sured in ordinary experiments. Then 

pCpSff ~ SQ ~ pcSff + 

~ pchd -H ^kadhij^hcijt 
= pcSd + ZkaO . 3aSd. 

^ka’^9\ __ k'c 

> “v* ■ “tt; “t 

The equations of motion at constant temperature, if the 
properties A and are uniform, can be written 


Thus 


= c( 


1 + 


(56) 


pfi ~ ^ (ASffcA -f- 2p,eift) + pXi 

^ dXi^ ^ ^ dx^ \dXf ^ dxj ^ ^ ' 

*= {(A -h ;i) A} -h + pXi 


(57) 


(58). 

JCT 
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If the changes are adiabatic A must be replaced by 
A'. 

If there is any heat conduction, the absorption of heat 
per unit time by the element of volume dr is 



where K is the thermal conductivity. Then the equation 
of heat conduction is 


Zd TZnoQ 

P‘^di + ~dt = ^ 

(59). 

If we write 

+ e,/ 

(60), 

so that e,fc' is the strain due to the stresses, 


dt dt 

(61). 

and the equation becomes 



(62). 

The element dr originally had volume 

0 

0 - Uj , ara - Uzy ara - Wa) 

0 (Xj, ^2, a^s) 


-{‘-(S+S+S)*" ''">■} 

dr (63), 


so that its density was p (1 + A). Hence 

dt p dt 


This is the equation of continuity. 



CHAPTER IX 


HYDRODYNAMICS 

In comparison with a typical elastic solid, a real fluid shows 
a great resemblance and a fundamental difference. The 
force per unit area across an element of surface parallel to 
a coordinate plane again constitutes a symmetrical tensor 
of order 2, and for the same reasons. The equations of 
motion are still (9) of Chapter vin, and also the rate of 
performance of work on an element of volume dr still has 
the form (16). The difference is that the internal energy in 
a fluid does not depend directly on how much it has been 
deformed. However the fluid is moved about and stirred 
up, provided it returns to its initial position, density, and 
temperatiire, the initial and final internal energies are 
equal. The deformation, however great it may be, makes 
no contribution ; the stresses do work on each element, and 
thereby supply energy, but this is removed in restoring the 
original temperature. If energy of deformation existed in 
a fluid, all particles of it would have a tendency to return 
spontaneously to their original positions when stresses are 
removed, and they have none. Accordingly, while the rate 
of performance of work on the element of volume dr in time 
dt is still pjciduijdxj,dr, where the symbols have the same 
meanings as inelasticity, we can no longer assert from this the 
equation (18) of Chapter vm, because the change of internal 
energy is not determinate when the changes of the are 
given. The may be as great as we hke, but the energy 
does not increase indefinitely apart from changes in density 
and temperature ; and the fluid moves in the same way under 
the same external forces whatever its previous history. We 
may say that an elastic sohd has a memory ; a fluid has none. 
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The stresses in a fluid are related, not to the total de- 
formation, but to the rate of increase of deformation. In 
the former notation these have components 


0e.- 


ik 


or 




du( dit 


+ 


dxj' 


dt ^ ' dx, 

The velocities now appear, instead of the displacements 
from an initial configuration, and we now denote the 
velocities, instead of the displacements, by 
We also write 


^ik 




du 

VBa:* ' dx^ 

so that now denotes the rate of increase of strain and 
the local angular velocity, usually called the vorticity. 
Now we say that p,* is linearly related to the new e,*; and 
therefore 


Pik — "h mp®mp (2)* 

where is a symmetrical tensor of order 2 that possibly 
does not vanish with the e,* , and is a tensor of order 4. 
Further, in a fluid at rest the stress is an isotropic tensor. 
Hence is isotropic and must be of the form — 
where p is a scalar, and must be of the form (20) 

of Chapter vii. On carrying out the summation with regard 
to m and p the term in v disappears and we have 

= — pS,-fc + ASjfcemm + (^)» 

where A and p are scalars, at present undetermined. 
Evidently p and A are not both necessary to preserve the 
form, and we may introduce the further convention 


Pa = - 3p (4), 

so that — p is the mean of the three with equal suffixes. 
This gives 

— 3p = — 3p + 3Ae„„, -h 2/xe,( (5), 

and therefore A = — (^)* 

If we consider the internal energy as a function of the 
density p and temperature d alone, we may consider the 



HYDRODYNAMICS 


85 


energy change in a symmetrical expansion. If the density 
increases by there is a contraction in aU dimensions in 
the ratio i^p/p and the stresses do work 

(’)• 

At the same time there may an absorption or generation 
of heat ; then the energy change is 


SE = p^ + SQ 
P 


( 8 ), 


and the heat absorbed may be written 

SQ = MBp + pcSe (9), 

where M is unknown and c is the specific heat at constant 

volume. Then x 

BE MjBp + peso (10), 


8Q jy pc 
-0 = + -0 


( 11 ). 


and the condition that these quantities are perfect differen- 
tials gives 

‘ (9 


M= -e 


dp = 


dd\p 


(13). 


de^Kp 

It appears also that p must be a function of p and 9 
alone, for a given material; it does not involve dpidt* We 
may call p the pressure. 

* This amounts to saying that there is no dissipation of energy in a 
symmetrical compression or expansion. This is true in a gas on the older 
kinetic theory; hut Enskog has shown Svtnska Akad. Hundlingarf 

63, no. 4, 1922, p. 18) that p can with greater accuracy be given by 

if Q.V^P 

where is a ** second** coefficient of viscosity. But tiIh is only of the order 
of the square of the ratio of the volume of the molecules themselves to the 
whole volume of the gas. 

In a liquid the coefficient of dpjdt, if it exists, is within the experimental 
error; an analogous statement is true for imperfectly elastic solids. 
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The stress components may now be written 

Vik = - (i> + f/xemm) ^ik + 2/lC(fc (14), 

where 2 ? is a function of the density and temperature and 
fjL expresses a property of the fluid. We call y. the coefficient 
of viscosity. For a unif orm fluid we have the equations of 
motion 

dUi , du 


= A 1“ (fe + 


dXi )] 

0 0 


dXi 


■= - 0— (v - 

where we now write the scalar 


(15), 


'mm 


- ^ 


(16). 


In time dt the outflow from a volume element dx^ dx^ dx^ 
is {d(pUi)jdXi} d/rdt. The mass within the element therefore 
decreases at a rate {3 and we have the equation 

of continuity 

dp ^ 

~dt “ dXi 




(17), 


or 


dp dUi . 

Ir-Pdxr-p^^ 


(18). 


Here r denotes differentiation with regard to the time 
dt 

dx 

following a given particle of the fluid, so that are 
given by 


and 


dXi 

-di = 

dt dt^ dt dxt dt *dxt 


(19), 


( 20 ). 


In most works on hydrodynamics this operator is denoted 
by DjDt, but I see no reason for departing from the usual 
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notation for total differentials, since this particular one 
is the only total differential that occurs in hydrodynamics. 
The acceleration components are given similarly by 



Consider now the circulation Q about a closed contour. 
We have 


S - s - IJs + L*'*' 

-L(^' 


+ - ) dXi + f UidUf 

P OXm/ Jc 


( 22 ). 


The integral always vanishes; X,- dr,- vanishes 

Jc Jc 

for bodily forces derived from a potential. Then 

= - d ip + f^A) + j J (- + 2 Ig)} rfx. 

= _ IJ d (p + j^A) - JJ gl; (Pl..) dx, 


+ f -p^du,+ f (23). 

Jcpdx^ )cp 

If jtx = 0 this reduces to the circulation theorem of classical 
hydrodynamics. In many problems of real fluids /t/. is 
small and constant, and A small. Then the integrals 

f -d(faA), f -|^di4fc,andf ? A are zero or products 

JCP JCP JCP 

of two small quantities, and can be ignored. -dp is zero 

JcP 

if p is a function of p, which is true in many cases, though 
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exceptions arise when the temperature or the composition 
varies from place to place. In the commonest case, to 
considerable accuracy, 

* ^ ~ Ic ^ 

and circulation arises only from variations of vorticity in 
the neighbourhood of the contour. It follows that in a 
fluid originally at rest or in irrotational motion circulation 
can arise only through the diffusion of vorticity inwards 
from a boundary. 

The work done on the element dr per unit time exceeds 


the rate of increase of kinetic energy by 



dUi , 

Pki 0 ^^ dr = Pik^adr 

( 25 ), 

and 

Vik^ik = — (i^ + + 2^e,fce,* 



^ - (P+ §/^^) ^ 

( 26 ). 

But 

1 

> 

11 

( 27 ), 


and if p is a function of p (including a constant as a par- 
ticular case) this is a differential with regard to the time 
and yields on the whole zero if the original density is ever 
recovered. This term may therefore be considered to give 
the increase of internal energy due to compression. The 
remainder may be written 

0) = - f^A2 + 2/xe,.fcC,t 

= {(^22 - 633)^ + + (Cii - 622)2} 

+ 4 /x (6232 + 63,2 + e,Y) ( 28 ). 

Thus O is analogous in form to the distortional strain 
energy of elasticity (( 47 ) of Chapter vm). It is essentially 
positive, and therefore represents work done on the fluid and 
continually stored as internal energy. In (8), therefore, if the 
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change is an actual one, the term (27) contributes thepBp/pj 
while (28) contributes to 8Q. So long as the initial and final 
states are given, it is immaterial whether 8Q represents 
heat conducted into the element, absorbed from radiation, 
or generated chemically within it, or mechanical energy 
dissipated into heat by viscosity. 

We notice that O can vanish if, and only if, 


ej 3 = 63 , = eu = 0 (29), 

SO that the deformation represents a symmetrical expres- 
sion or contraction. 

If we consider any finite volume, the rate of dissipation 
within it is 


But 

and 


in = 111 (- + 2 iiei^eit)dT 


, ^ 


(30) . 

(31) , 



where is a direction cosine of the outward normal to the 
boundary. But 

= II 2l„Uii,kdS + i 11^’ 

= ||2.,,„«Af„d5 + i||f-'d5 (33), 

where d/dn denotes differentiation along the outward 
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normal, q is the resultant velocity, and is the vorticity 
vector associated with f,*. Also 


J 4 


u 


*dx 


dxf 




and in all 



Odr = ^ ^ + 4/x j| €ij,„UilkiyndS 

+ + + (35) 


This form is useful in such a problem as that of waves 
on deep water, where the viscosity is small; if it were 
absent we should have a permanent oscillation in a normal 
mode. The vorticity is negligible except within a distance 
from the bottom of order (v/y)^, where ft = vp and 2njy is 
the period of the motion. Then in the second integral 
is zero at the free surface and Ui is zero at the bottom, so 
that this integral vanishes. A is negligible everywhere. 
Near the bottom contains a factor proportional to the 
velocity that would exist there on the classical theory, and 
thus can be made indefinitely small for deep water. Hence 
the important term in (35) is 

<• IJ S «'»■ 

and even this vanishes at the bottom, so that it need only 
be estimated at the free surface. 

If a portion of the fluid is compressed without change of 
temperature, we define the bulk-modulus k by 


1 

p dp k 


(36), 


and if it expands under change of temperature without 
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de~ " 


(37). 


change of pressure we define a coefficient ol volume 
expansion a by 

1 dp 
P 

Then for small changes of temperature and pressure 

- 

Also in a free expansion under constant pressure 

= Mdp + pcdQ 

d 


8p = P (' 


(38). 


r (p) ■ 


ap + pc 




(39). 


where the partial differentiation is to be carried out at 
constant density. But this makes 


1 -“ 


(40), 


(41), 


{ k ) 

and therefore 8Q — p |c + ^ 

so that Cp, the specific heat at constant pressure, is given by 

a^kO 


c« c -I- 


(42). 


This form is analogous to that of (56) of Chapter vin, the 
present a being the coefficient of volume expansion, and 
the previous one that of linear expansion. 

In an adiabatic change dQ = 0; then 


— akO — 4- pcd9 == 0 

P 

and therefore aOdp = {pc + a^kB) d9 

a^k9\ ak'B 


(43), 


( 


= 1 + 


V V ^ 


pc / p 


+ “-f) 


80 that 


(44). 
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Thus the bulk-modulus for adiabatic changes is 

‘ + "^) - t’ !«>■ 

The equation of heat conduction needs to be modified 
to allow for heat generated internally. In time Bt, per 
volume dr, we have 

se - * |(ii: I) 8. 

= pcB9 4- akOABt (46), 

whence pc | + ak9A = <D + A (z g) (47). 


EXAMPLES 

1. Obtain the equations of motion in terms of the stress 
components by considering the momentum of a finite volume of 
any form and applying Green’s theorem. 

2. Similarly obtain equation (25). 
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Angular momentum, 36 
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Curl, 14, 68 
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line of action of, 42 

Gauss’s theorem, 49 
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Green’s theorem, 48 

Hamilton’s principle, 28 
Heat conduction, 82 
Hydrodynamics, 83 
Hydrostatics and classical hydro- 
dynamics, 56 


Parallel, 17 
Potential, 49 

Rigid bodies, 29 
Rotating axes, 40 
Rotation, 21, 30 


Scalar product, 7 
Scalars, 3 

Specific heat, 81, 91 
Stokes’s theorem, 48 
Strain tensor, 72 
Stress tensor, 73 
Suffixes, dummy, 6 
Summation convention, 3 

Tensors, 2 

antisymmetrical, 5, 8 
isotropic, 66 
symmetrical, 6 

Vector product, 7, 14 
Vectors, 2 

with given divergence and curl, 
58 

Velocity, 24 
Virtual work, 34 
Viscosity, 86 
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Inertia, moments and products of, 37 Work function, 26, CO 

Isotropic solids, 78 "Work, virtual, 34 

Wrench, 45 
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